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COLLECTIVE OSCILLATIONS IN AN ELECTRON GAS 
OF METALLIC DENSITY. THE CORRELATION ENERGY. I 


By JAN ŚLEDZIK 


Institute of Physies of the Polish Academy of Sciences, Department of Ferromagnetics, 
Poznań. 


(Received June, 27, 1960) 


The role of the subsidiary conditions in the Bohm and Pines theory of an electron 
gas is analysed anew. It is pointed out that the subsidiary conditions must be formulated 
in the form of two sets of operator identities: The first set establishes the Bohm-Pines 
operators of subsidiary conditions, and the second — the commutators of the extended 
Hamiltonian and of the first set. The extended Hamiltonian and the two sets of supple- 
mentary conditions are proved to be equivalent to the /V electron Schroedinger equation. 
The two sets of subsidiary conditions enable, at last theoretically, to replace a certain 
number of electronic positions and momenta by coordinates and momenta of a longitudi- 
nal field describing the long-range density fluctuations. This problem cannot be solved 
in practice on account of mathematical difficulties. Therefore the assumption is made 
that a certain number of excited electron-hole pairs and the collective response to 
them — the rearrangement of the charge — may be described by a certain effective longi- 
tudinal field as introduced by Bohm and Pines. In the theory thus formulated, the total 
number of degrees of freedom of the system is correct and no subsidiary conditions are 
are required. The Hamiltonian and the ground state function are the same as in the Bohm- 
Pines theory, but the number of individual electrons is smaller. The correlation energy 
obtained is several times smaller than that of Bohm and Pines. The theory proposed is 
valid within the alkali metal density range. 


1. Introduction 


In the present paper an attempt is made to elucidate what form the subsidiary 
conditions in the Bohm and Pines theory assume and what is their role in the collective 
description. This problem was discussed in several papers [1], [2], [6, 1957b], [12, 
1958a], but the conclusions were controversial. Bohm and Pines formulate the subsid- 
iary conditions in wave equation form [1], or in that of operator identities [3]. The 
former fails to account for the reduction in the number of degrees of freedom (in 
the extended Hamiltonian a certain number of new longitudinal degrees of freedom 
are introduced) to the value 3N, as we suppose the physical system to consist of N 
electrons. Therefore the field of the longitudinal oscillations does not describe the 


(3) 
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long-range charge density fluctuations, but represents some external fictitious field 
with which the electrons are very strongly coupled. The state vectors of such a system 
are non-normalizable, hence the perturbation method developed by Bohm and Pines 
in [1] is divergent. All this was proved in the paper [12, 1958a]. 

As we believe, only the subsidiary conditions in the form of operator identities 
are legitimate. This problem is investigated in detail in the following section. 

There are two sets of subsidiary conditions: the one relates the momenta of the 
oscillator field to the long-range charge density fluctuations, whereas the other one 
relates the oscillator field coordinates to certain functions of the electron positions 
and momenta. The second set is not independent of the first. The simultaneous solu- 
tion of both sets of subsidiary equations, i.e. the replacement of a certain number 
of electronic positions and momenta by oscillatory coordinates and momenta is too 
difficult mathematicaly, but this nevertheless is the only procedure ensuring that 
the total number of degrees of freedom of the system is conserved and simultaneously 
the long-range density fluctuations are represented correctly. 

Since the aforementioned problem cannot be solved strictly, we formulate a new 
approach in section 3, consisting in the replacement of a certain number of excited 
electron-hole pairs and the changes in density of charge they cause by an effective 
oscillator field describing the long-range density fluctuations. Since now the total 
number of degrees of freedom is correct, the subsidiary conditions become superfluous. 

A beautiful picture of collective oscillations expressed in terms of excited electron- 
-hole pairs and individual electrons is due to Sawada, Brueckner,, Brout and Fukuda 
[11]; their theory, however, is valid in the high density region only, whereas that 
presented here refers to the alkali metals. The constant in the correlation energy 
derived is two times smaller than that obtained by Gell-Mann and Brueckner [5] 
and can be further improved by perturbation treatment, since we used the ground 
wave function of uncoupled plasmons and individual electrons. 


The subsidiary conditions in the Bohm and Pines theory 


It may be useful to recapitulate the procedure adopted by Bohm and Pines in 
deriving their theory [1]. Let us assume a system consisting of N electrons. Instead 
of introducing the positive ions, which complicate mathematicaly all the computa- 
tions, we assume the electrons to move within the volume V filled out with a uniform 
positive charge background. The total charge of the system is zero. All the calcula- 


tions in the present paper are valid when N > + oo, V = I3 > + oo, — = const. 


V 


We assume the volume V to be a cube of edge L, and apply the periodic boundary 
conditions. 


'The Hamiltonian of the system is now 


x N 
f ; I Ne? dr N?e? dr d 
GIER p E 4 Ne e rdr” 
2 DESEE e tire y X ZE 272 fe P|’ (21) 
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wherein the third term represents the interaction energy of all electrons with the posi- 
tive charge background, and the fourth — the interaction of the bac kground with 
itself. 

The Fourier transforms of 1/r (r is the position radius) and Dirac's function are 
of the well known form 


l 4x elk: r 
pm ap CZ 
k 
1 s 
8()— D enr. (2.3) 
k 


It is also convenient to introduce the Kroenecker delta, 


k=, 


0, otherwise. 


à = | (2.4) 


The charge density of the system is now 


N 
Ne 1 ! 
e(r) =— eY ôlr r) + y Sy, eet, (2.5) 
zl k 


where Ne/V denotes the density of positive charge background. With respect to (2.3) 
and (2.5), we obtain 


N 
— e Y etti + Ne 6, (k), (2.6) 
i=1 
since 
Ne Ne HR 
Wy 2: dx (k) eit n. (2.7) 
Moreover, 
N : 
09 = — e 9] + Ne = — Ne 4- Ne — 0. (2.8) 
i=1 


Let us now transform the terms of the Hamiltonian (2.1) one by one. For this 
purpose we use (2.2) and vod 


> 1 eż x ce etk: rjj Sr S "Y 
i ates m Lm 


(i) lij iz i£j k i,j 
la N e? eV 1 
y ji 2 (2.9a) 
k 
on the other hand, 
N NU 
0,0, = [— e >; e 7 + Ne dg (K)] . [- e 2 e’k'ri + Ne dx (K)] — 
Um i=1 
=e? DE "=. Ne? 0c (Ip), (2.9b) 


5j 


or 


2z po ko xum cik’ rij 2r N? e? Ók (k) 
V - D a [7 a k2 ? (2.9c) 


hx k 
By comparing (2.93) and (2.9c) we obtain 


' | 24 Ok0-x — Ne? , 2a N? e? Ox (k) 
ee) etme ja ey > EA (2.94) 


Te 
i<j K k 


The third term of eq. (2.1) is now 


327 — br i E eik- r 


i=2 * i 


Ne 4a N? e x Ox (k ) 
= fes a SYL - V dx (k) = — > aus (2.9e) 


k 


Quite similarly, the fourth term of (2.1) is transformed as follows: 
2e ('l' dr dr? 27 = e? b 
2 E ; > Y = dr eik'r | dr! eir 
ay? JJ r=r| k , 


2a N? e? | PZA 27 N? e? Ox (k) 
LL - m r >») — ——-_. _— E 2 
V3 » ja 5x (h) V M : en 


k2 
k 


Collecting all transformed terms we get 


N 
A Di 9 = 97 N? e? N? e? 
=) y EN: _p—Ne? a S 


1=1 
„że ye sA NET Seen Nels ie cae 
V De: Qin V = ENT j » NS 


This is modified form of Schroedinger's equation for the system of N electrons and 
positive charge background. 


Such a system of particles may also be described by the longitudinal field repre- 
sented by the vector potential 


dar c? * ode > k 
Awa ^y- Da m (2.11) 
|k| Sk, 
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where c is the velocity of light and gq, are coordinates of the field. The electric vector 
is now 


= Wr Ar Nc 0%. | 4 
E(r)---——4G = Jk 5, eik'r — 7) d, eik-r 
W SER ZAŁ V pia -y5 TEE 


|k|<k, [k| < ky (2 12) 
with p, as momentum of the longitudinal field. 
With respect to (2.11) and (2.12), we have > 
Ê l x A e Er 2z Ok O k — Ne? i 
|" 24 ALS a ifo mA = Ha (a 
uem n Dus ds) c jà peeo- 
ta |k|>k, 
AIN NT ol 
E 3 a (2.13) 
|k| & ko 


We postulate now that Maxwell's equation be satisfied: 


[div E (r) — 4xo (r)], = 0, (2.142) 
OT 
D= 1/7 o = 0, |k] < ke. (2.14b) 


The subsidiary conditions (2.14b) relate the momenta p. , to the long-range Fourier 
transforms of charge density, i. e. to the functions of all electronic positions. Let us 
write the conditions in explicite form 


RN 
4 
QP = P-k + 2 ES da Qr S Os alk s. (2.15) 
i—1 


Were we able to solve these equations, it would be possible to replace the n’/3 electronic 
positions (m is the number of all k within the interval |k| <k,) by m oscillator field 
momenta, but such a solution of the problem cannot be obtained. 

Let us now rewrite the Hamiltonian (2.13) in the more developed form using 


2.11) and (2.12), 


N 


^ Dicen Dat om Ne 1 
Ss SS ps E > e n (pa + 02ga gx) + 
i=1 Ikl>k, ~ |k|<k, 
==), 
e As À m 2a Ne? 1 
FE zV DE qr Er" (Pi — hk/2) en = pg 
i=1 |k|<k, |k|<k 
N 
2 = 4a N e? t 
ad 27 e Da qe die te e V | elk Tk’) r; w? — id: y 5 (2.16) 
DE ise, a is 
|k?|<k 


It must be 


99] = 0. 


(2.17) 


[H new ? 


This is legitimate, since QV, according to (2.15), is the zero operator. Straightforward 


computation yields 


di Qa = 2ieh m ŻAL (f; —hk/2) eit vi 
Ar 
er | r dy Ex‘ & eet | BU (2.18a) 
or 
. N E 
OH Ep (Di — hk/2) eit v 
ri N 
T wt 
T VF qe Er" Er eet) = 0, |k| e E, (2.18b) 


This is the second set of subsidiary equations relating the oscillator field coordinates 
to the electronic momenta and positions. The solution of these equations, if achievable, 
would replace the n'/3 electron momenta by the longitudinal coordinates q, (we assume 
the n° selected r; to have been replaced by the momenta p,). In this way two sets of 
subsidiary conditions would yield the 3N — m electronic and n° plasmon degrees 
of freedom. But this programme is mathematicaly too difficult for execution. 

The physical equivalence of the Hamiltonian (2.16) together with the subsidiary 
equations (2.15) and (2.18b) and of the Hamiltonian (2.10) is easily proved. The proof 
is along the following lines: 


Ag 1 An 


R=il/ pe- Pa=i]/ pow, l|k&k,. 


whence 


l 27 A | 
SUED Y PA Pei p M AS. (2.192) 


By (2.18b) we have 


N 
Du 2 * (Pi — hk[2) een; + e in 
x] 


N 
dk Èk * Eh N ee N 


|k’|<k, i=1 


Multiplication of the last equation by yF q, and summation over |k|<k 
m 
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(the double sum over k and k’ must be multiplied by 1/2 to avoid repetitions in K) 
yields 


== IN! 
e 4n AA E e? 
>V V dk Ek * (D; — hk/2) ek vj + >) dk Qi Ep: Ep > gie k?) * r; 
i=1 |k|<k, 


Mers i=1 
|k’|<k, 
—N 
e D X (p hk|2)eik Y 
c qe = dk Èk ` qim e ta CO dk 
m V A ee, 2 RSR, „p 
N 
ONU EE > RR — 0, (2.19b) 
|k|<k, = 
L'ART W 
k+k' £0 


With respect to (2.19a) and (2.19b), the Hamiltonian (2.16) is of course identical 
with one of (2.10). 

In the paper [1] Bohm and Pines postulate the subsidiary conditions in form of 
the wave equation 


QUO IS he (2.20) 


where Q is given by (2.15). Such equations are satisfied by non-normalizable state 
vectors, as proved by Kanazawa [6, 1957b]. Moreover, the principal requirement, 
which consists in the reduction of the 3N + n’ degrees of freedom to 3N is not ful- 
filled by such subsidiary conditions, hence the solutions of (2.20) describe the be- 
haviour of N electrons in some fictitious external field. The interaction of electrons 
with such a field is very strong, as computed by the present author [12, 1958a]. Thus, 
obviously, the perturbation-theoretic treatment of the problem given by Bohm and 
Pines [1] is fruitless. If we require only 


QD —0, |k|«k, (2.15) 


the theory is freed of aforementioned difficulties, since now some electronic degrees 
of freedom are replaced by other ones, which also are internal. The formulation of 
subsidiary conditions in the form of (2.15) was investigated by Bohm [3]. 


3. Collective oscillations in an electron gas. The correlation energy 


Suppose we have N electrons immersed in a medium of uniform positive charge 
within a cube of volume V. Let us at first assume that the electrons do not interact 
and that they are regularly distributed over the entire space qf the system. By regu- 
lar distribution we understand that any one of the N electrons is nearly at centre 
of its” cell VIN. The Fourier charge density transform is now nearly zero. Let us 
now very slowly (adiabatically) to *switch on" the interactions between the electrons 
until they attain their real values. Since there are locally slight differences of charge 
density, these increase with the increase in interaction. Finally a certain number of 
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electrons enter into collision, giving rise to local charge unbalance imitating an excited 
electron-hole pair (the “hole” is here the unbalanced locally part of the positive 
charge background). The system as a whole meets the appearance of such pairs with 
collective response, i.e. with some kind of rearrangement of charge. Both interconnected 
phenomena, namely, the arising of a certain number, say n’/3, of pairs and the new 
distribution of charge depending thereon will be described with the help of an effective 
longitudinal field given by the vector potential 


Air) ae > > qy & eif r, (3.1) 
HER, 


We assume that the state of the system is such that there are N — n’/3 individual 


electrons and n’ field coordinates q, responsible for the n'/3 excited pairs and for the 
related rearrangement of charge. We further assume the state of the system to be 
“stationary”, i.e. that n’ depends only on the temperature and the mean distance between 
the electrons. 

We now write the Hamiltonian of our system: 


N’ 
^ il A ee 27 Ok 0-4 — N'e? l z 
== ; 4.) 208 = = ~ Ffy) — 
z 2m > (p; + c A) V 3 k2 8a fo BATA 


1-1 |k|>k, 
_2aNe 1 Ne S | dt SO Ne] arr ah 
V BUY r—r wejJr=r" 69 
I| Zk, i=1 * | 
wherein 
> 4 | 
p E (1) = V* M pré e^t (3.3) 
|k| ZR. 
Moreover, 
N’ = N — n8, (3.4) 
» = >). (3.5) 
|k|<k, 


Ihe second term in eq. (3.2) represents the short-range Coulomb interaction between 
the electrons, since it may be additionaly assumed that the rearrangement of charge 


due to the excited pairs produces screening in the Coulomb interactions of the indi- 
vidual electrons. 


The total number of degrees of freedom is now correct: 
3(N — n’/3) electronic plus n’ of the field = 3N. (3.6) 
In equation (3.1) we include k — 0, 


Pe M 4x c? pate oca Age c? 
A) s] PEL. M dk Er eh "y y D: (3.7) 


PETS 


where prime in the sum over k denotes the dropping of k = 0, and the constant qq 
will be determined later on. 2? is the unit vector along the z — axis, 
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By (3.2) and (3.3), the explicite form of the Hamiltonian (3.2) is now 
klr;-r;l 


N° ^9 
Bi 2 n 2 sin T 1 
"PEE 2m tey li? | a M 
i<j 0 |k|<k, 
rrp 
2 € TT s5 + > 5 
r Wp Ik q—x) T PX HE Ds D dk Ek’ (Di = hk/2)eik'r; + 
i=1|k|<k, 
2a e? = Ast e? aL , 
>> = = SA oh: BUBLE = $ 
+ my $ Y dk Qj Ep ` £e ETE) r; 4 my do Di » gp 20 + £j ek" ri + 
i=1 |k|<k, i=1 |k|<k, 
|| E, 
ktk’ £0 
2a N'e? 2a N'e? UNE 
2 a 
re = Se =) Ina Frj 
|k|< ke 
dr an 4g N’e? 
| 02 = 
| e eff a een Fe) 
We write the ground state function of the system as follows: 
Pr, ol), 2; Pa (nr) Ny (On) 
3 1 i Te (ri) ny (CAES Pry, (ry) N+ (Sn) 
wi : (3.9) 
ya 9-9) 1 (Gw) 
Pen (r1) n- (03); +; Pan (nr) N- (on) 


where ®, is the oscillatory vacuum vector, 7+(0;), 7-(0,) stand for right spin respectively 
left spin functions, n, is the number of right spins, and 


1 | 
qx (r) = ZA CUR (3.10) 


We shall investigate the two magnetic possibilities only, namely 


Ÿ = 1/,, non-magnetic state 

E Ie e (3.11) 
Ÿ = l, magnetic state. 

In connection with eq. (3.6) we introduce the Fermi distribution for absolute 


zero temperature 


y. A = (27)? N°, (3.12) 


wherein kę is the wave vector of the Fermi surface, 
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The Fermi distribution for the N electron problem (without plasmons) is of | 
the form 


Ast 


V- kb (2m)? NO. (3.13) 


ką now depends on N. 
The number of longitudinal degrees of freedom n’ will be 


|k|<k, |k|<k 


ke 
y V RAM 
P = = ap LA paul | an Rt a sE 
t 2 Gas J 4 eal St ae E 
0 


The quantity n’ or k, is the variational parameter of our problem. It will be more 
convenient, however, to solve the variational problem with 


p= teu, (3.15) 


and not directly with k,. 


By (3.4), (3.12), (3.13) and (3.14), (3.15), we obtain 


AR n Maj mp 8 
Sw = U ZR 

N : 3N | N 2:31 E 3 p, (3.16) 
kr Ü p LA i 
ko ( TE 2 (0047) 

We now introduce the spacing parameter r, by the equation 

47 : l 

EA (3.18) 


with a, denoting the Bohr radius, and by combining eq. (3.18) with (3.13) we get 


L Jad |*1 me? 
Uu) TE 2 T, h2 . (3.19) 


Let us now compute the mean kinetic energy per one individual particle in state 
(3.9) and on the average over the Fermi distribution (3.12). We have 


1 N ng h2 "ża 9 42,2 
T — Ne p pi " NG (of [dk k 3 h? kp 
N’ , — 2m d 2m [dk D Se Du 


_ 3 h?* o (kp)? 3 [920 " 1 Dam % 1 met 
5. 2m ES] DS er r: 2h? 


3 [920^ Duc EY 
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We have utilized (3.17) and (3.19) and expressed the energy in Rydbergs. For the 
non-magnetic and magnetic state we obtain, respectively, 


2.2100 1 % 
Tes 2 | 1 p) Ry, (3.21) 
ye 6 
3. 5077 s 
T > SE EJ bat ka ie) Ry. (3.22) 


In order to compute the mean of the oscillator part of the Hamiltonian (3.8) 
we notice that 


he, |^ : 
Der s (a-r + ak), 


a 


h Vs a 
J= ER (ar cm GE k = 0, (3.23) 


* . . . . . . 
where a,, a, are annihilation and creation oscillator operators. Since the vacuum 
vector is defined by equations 


GoD KEN (3.24) 

we have 
DT D 200. ID, p. dos 0, (3.25) 

Moreover 
Hose = — 3 D ip- + esae i u 2 (8:20) 

whence 
1 A 5 ho, ho, , 

ip s Do(PrP-r + pqr 1-4) Po = » E Pur (3.27) 


|k|<k, Ik|<k, 


By (3.14), (3.15), (3.16) and (3.17) we obtain per one electron, 


1 p ZA ^ ho, n 
ivy wO Ales v = ys Pó Hose Po 2N° 
Va he, (N\ hoy 8 ga 
N 2-3!g81N' 2 sł 

1-28 


E. expressed in Rydbergs is 


osc 


Boe == 


(3.28) 
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The quotient ho, ag/e* can be transformed as follows: 


ho,ayM  h?4a Ne?żag _ EB! T e iu 1 — Rare 
(= | EN ag a LES ne qu 
— agr; 
3 
3 Ding 
Finally, 
ho, V3 1 U p) 
e?|ag — rh 3 
and 
3 
E, de SĘ R 
osc 3/ Ü 1 
s a= 3 
z 
Thus, we have 
3 
di = £ P y, Ry, 3.29 
aris 1 Lgs| * (3.29) 
6 
R 3 p? 
att z Ry 
ls RE 1 (3.30) 
3 


The mean value of the fourth, fifth and sixth terms of the Hamiltonian (3.8) is 


zero owing to eq. (3.25). The zero component of the fourth term 


N* 

PEN da DI 29. p. 

m y 4 í P: 
i—1 


vanishes since the Fermi distribution is isotropic in momentum directions. 


> 2 
The long-range self energy term — zy 1/k? is, per one particle, 
I| A, 
w 2n e? INS e? dk a NE e? ko B 
y V k Agit] A It T 
|k|<k, Ik| Sk, 
«2 il B met 2 [9x0 |^ B 
mE, FORME TE NI d 

L2217 

ET == EU Ry (3.31) 
1.5392 

E) =— p Ry. (3.32) 
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Let us now examine the screened Coulomb interaction term in (3.8). For this 
purpose, we compute the following mean value: 


k,|m=rr| 


2 s] 
(w. € ( 2 fa a v). 
In, — rej 7 T 
0 


By (3.9), we have, 
1 c Pk, (r3). Pk, i6 2) 
k p vi, (Pa), (r) 


> Pr, (11), Pr (2) 


Pr, (rj). Pky (r2)| 7 
CNE) | 

| 

| 


t 


wherein the 312 denote the minors (including sign) of the determinant (3.9) obtained 


by deleting therein the /-th and p-th lines and the first and second columns. Summa- 
(GM) EE ED) 
tion >} or >} runs over right or left spin indices. 


On account of the orthogonality of the functions (3.10), we obtain 


(ps 12) 3 = 2 di dr, kai dr w yp*lp;12 yuv; 12 


My. JU 


= (yu ctor 


4 (01) 1. (02) PR + 


1- (01) q— (0) FB: 12 + 


Pk, (13) 9. (03), Pk, (r5) N4 (0 


Pky (1) N- (01); Pe, (Ta) N- a Pon), 


dT dep, cu ŚW, (3.332) 


hence 
km: 


e? 2 sin T 
* ee y 
E: nr; 73] E A fa T | | 
ô 
(+) 
LX jh Ert Ji ins Y 


vi (rs), Pi, (2) 
Pky (r1). Pi, (2) 


Phy (ri), Pky (rg) 
AE (11); Pk, (r;)| * 


ne) m (03) gn is atc 


x mh (01) 13. (0) bs + > 


ke|vi — ra| 
á é 
ep; 12 e 1 2 E sin J © 
[ry — rz] d Ę 
0 


Pky i Pr, (12). 
> Pr, (T2) 


nn Pk, (03) n+ (01) Pk, (ro) 1% (02) 


gk (11) n= (04). qx, (ro) N = (02) 


Pr, (11): Pr, (T2) 


Pr, (ri); Pr, (r,)| "+ (03) 1, (02) FP: 12 + se 


jo (01). 1- (03) X 


Jb» 
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CLS 


x Wuv; 12 PRZ 
x Sx 102 p Jar 


I<p uco 


Pk x r1) N4 (01). Pr, (72) N4 (03) 
A 1)» 9, (r2) N- (02) 


Wuv; 12 — ea (SS 


I<p u<v 
Rafis 


= |= 
12 

equ sin T GE 
EL IE > N, p far f ar, toż, o (3) Pr, (T1) Pk, (r2) Pe, (Fa) + 


| Pky, (1); Pry (la) 
Pr, (1). Pk, (r3) 


Ph, di vi, (2) 2 
> Pk, (ry) 


0 
Rove 
* * e? 2 sin T 
uis Pky (rj) Pk, (rı) Pk, (ro) Pk, (r3)] Lr i= = i dt E 
0 
k.lrı—ral 
dr, dr, 2 sin T 
T m 
cei oS ARES a [^ | 
0 
(+) (-) k|r—rl 
e? egi - hj) *(rı— ra) 7 sin T 
ma zel ‘2 3335 
> px T Y) [far ar. AEST ( = | * - ) (3.33b) 
l<p l<p ` 


In calculating (3.33b), we made use of the well-known fact that the spin functions 
are orlonormal, 


M he (c) n+ (o) = > n (m) »- (a) 5m 


D Na (0) n-(o) = 2) »- (6) n+(o) = 0. 
Since M- MESS mj , we have, finally, 
l<p 
kom 1] k.ir—r’| 
2 i b ( ed 2 L E |= AŻ dr = erdt n sin T = 
[ri—13| Ite. T r- TT t 
0 0 
klr=r'| 
> | ws eilky—kj) (r-r) 2 ex 
a oe dr dr? í ; l [or (3.33c) 
pr — 2| a T 
l<p l<p « 
0 
Or, 
ker 
DIE: e? 2 : 
T y= l_— — [dr sinv| V 
u jj T 
i<j T 


0 
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k.lr=r'| 
DELE "dr dr? 2 d sin T 
2y? Fr] ee) 27 
0 
k|r—r?| 
GE) c 
is =k) (rr) 2 sin T 
vs 103 n Y) IET: 7 ( - L : ). (3.33d) 
I<p I<p 5 


In the last term of eq. (3.33d) we make the substitution 
1 > 

oe iS dr dr” = dg dR, 
2 is 


> 


"e N 
and proceed to the limit N’ — + oo, V — + 00, 7 == const, 


k.lr—r'| 
eil — kj): (r—r") 2 sin T 
a? uee 
ff eo Fr] ( JE x | 
0 
k 0 
i(le, k) © sets: 
emn 2 | ac) fan 
0 v y 
0 
+00 Foe 
Ag V 4 2 sin T 
a PT = 1—— [d 
0 0 
k,@ 
AV 1 | 2 id Vi 
= silk, — k 1—— | dt — = 
xu ik, — Fy cos |k, i| o z. T 0 


ES qug e Any h 1 Y 
( 


o k, = - k)? TT 
T: Anm V 
T 
z UN (ke + |, — kil) o + sin (ke — |k, = Hp al} = LL see dl es 
3 (Kp 1) 
; (3.336) 
If we now define gg from (3.7) as 
mes ko? , 
2z N’ e? N (N' — 1) e? (far dr’ sa NON ye? (ide sin 7 
qĄ = | ; dt = y — | dz sza 
7 ek zy? ko ZAJ T 7 217 
Te. 


ke 
, 23 s] 72 
xi [S LES fe) (3.34) 
0 In? o T 
: 0 
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then the first term in eq. (3.33d) cancels the last two terms in (3.8) and the expression 


s 2 . LI 
a de. The last elimination is due to (3.34). The proof is straightforward: 
m 


, NQ- De [ dr dr” jer Jab dro. Ne I dr dr” 
Mec ZE) Je T TOES r—r] 272 JJ Fer 


WADA Lok Rare [E dà 
ET sera brin le 


Per one particle, 


l e? |do 
wt = =, ; > 0 (3.35) 


as V > + oo, since V increases to infinity as Z? (L — the cube edge), and the integral 


[= [acc 
£ 


We have yet to average the exchange integral of (3.33e) over the Fermi distribu- 
tion obtaining, per one particle, 


(+) 
Ar e? 
Ban — sy » + >) < Sc —Es = (3.36) 
I<p I<p 
wherein 
o char Fan 6 
Sm: f dk, f dk; + [Ky — ki| > ke, ko < kp. (3.37) 


The mean value (3.37) was computed by Pines [9, 1953a] as 


pe 1 »-ixh-$ksl i TE 
(ki? 4 3 kg ' 2 \kp 48 | kp 
1 
2 


2 RUS 5 „W |. (3.38) 
wg 
It is readily verified that 
(+) (Ga) 
NAD 
) + = (3.39) 
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whence 


In the two magnetic cases 


0.9163 TR p? pa 
He = In ) == | SS | 
i ^ 6 imi zo veo Otel 


DERI D Nt. 4 js p* 
Eexch = > | 3 p) SZA 2(1—183*^  48(1 Pus] pu 
(3.41) 


Assembling all results obtained (3.21), (3.29), (3.31), (3.40), and their parallel 
spin analogues, we have 


Tate 0.8661 f3 1222178 
2 | — 
DT ( 22 HRS SS rs 


0.9163 i eum a 6: E 
"e [(-5e) gon 20 —183)4 sat Ry, (3.42) 


75 
3.5077 (1- 1 p) xiu 2321.02 1.5392 f 


CHOSE on 


s 4 


1.1544 ed B 
ae [(- ze) GA BE mus Le À 


The relations (3.42) and (3.43) must now be minimalized with respect to 5, 


hence 
9E. _ 0.9163 8 215982) Wp i 015268 1) 2 10.7367 | FE 
ade 85/58 cnr (Led) u ee ra)” | 


0.0764 83 0.0764. 8% | 0.2165 65 0.000589 —— 0, (3.44) 
ET Lost P9 re ERIT nip on 
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2.3089 
20-16) 


1.15448 


de 


„(1 — 1897 


0. EE p? 


3.4642 
de ien 
Ur ( 


1-488 


0.1924 f 


+ 5°) 


5 (1— 3 psy * 


0.8848 —— 
RASY 
0.8661 85  , 0.024185 
n*(1—189'^  nü—ig) 
(3.45) 


We compute the correlation energy for the non-magnetic state only, since the 
energy of the magnetic state lies higher and therefore this state is not realized. We 
define the correlation energy as 


ka = E, P E, 
2:91 0.9163 
p-2200 _ 
Ts Ts 


(3.46) 


(3.47) 


In Table I, the above listed quantities are compared for some alkali metals. 


Table I 
| o 
Metal Li Na K Rb Cs 
T ji 2 3.225 3.96* 4.87* 5.18* 5:51 
Bom 0.4450 | 0.5466 | 0.6380 | 0.6820 | 0.7305 | 0.7454 | 0.7625 
ua — — 0.6536 | 0.6753 | 0.6994 | 0.7065 | 0.7226 
Ey, Ry 1.2629 0.0710 | —0.0915 | —0.1092 | —0.1126 | —0.1117 | —0.1103 
Ry 1.2937 0.0943 | —0.0715 | —0.0905 | —0.0950 | —0.0945 | — 0.0933 
Ey Ry — — —0.0198 | —0.0689 | —0.0918 | —0.0951 | —0.0970 
Bm Ry —0.0308 | —0.0233 | —0.0200 | —0.0187 | —0.0176 | —0.0172 | —0.0170 
Hox eue E —  |-0240 | —0.269 | —0.308 | —0.319 | —0.333 
| 
The correlation energy is approximated to —13% accuracy by the following 
formula: 
Fore = — 0.0502 + 0.0622 In r, — 0.0132 :r,Ry, 3«&r,«6. (3.48) 
The corresponding quantity obtained by Pines [9, 1958c] is 
= — 0.140 + 0.0622 In r, — 0.0038: r, Ry for all r,, (3.49) 
while Gell-Mann and Brueckner (1957) give 
jor = — 0.096 + 0.0622 In r, + O(r)) Ry, 7, <1. (3.50) 


Here O(r,) includes all terms 7? and 7” lnr, from n = l onwards. 


Figs. 1 


and 2 illustrate the data of the Table I. 


In the region r, < 1 inclusive of the high density region r, «£1 the foregoing 
calculation does not apply. This may be seen from (3.44) and (3.45), where the coeffi- 


* The values marked with an asterisk have been taken Prom Pines’ 


paper (19534). 


21 


cient of f* points to anomalous behaviour just in this interval. The physical reason 
for this anomaly is clear; namely, we used the ground state function in the form of 
(3.9), whereas in the high density region this form in no longer adequate for describing 


/8 
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44 
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Fig. 2 


the system. Therefore it becomes necessary to apply perturbation technique either 
in order to improve the state function (3.9), or, and this is what the present author 
prefers to do, to eliminate in some approximation the interaction term in the Hamilto- 
nian (3.8). This latter procedure will be followed in a subsequent paper. 

The author thanks Professor Dr S. Szczeniowski for his numerous discussions 
regarding the Gell-Mann and Brueckner high density electron gas theory, as well 


as for reading and commenting upon the text. 
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In order to isolate the plasmons in the Hamiltonian derived in the first part of the 
present paper, the canonical transformation to independent plasmon modes is applied. 
The formula for the correlation energy obtained is nearly the same as that of the Gell- 
-Mann and Brueckner high density electron gas theory. The absolute value of the constant 
in the correlation energy is by 4 % larger than in the Gell-Mann and Brueckner approach. 

The second order energy shift within the Rayleigh-Schroedinger perturbation theory 
framework is evaluated. 


1. Introduction 


In the first part of the present paper we proposed a new formulation of the collec- 
tive description without subsidiary conditions. Thereafter we evaluated the correlation 
energy of the electron gas system. 

In the present part we proceed to transform the Hamiltonian (I 3.8). The aim 
of this procedure is to uncouple the individual electrons from the plasmons. Subse- 
quently we apply the standart method to obtain the ground state energy and the corre- 
lation energy. 

In order to obtain the second order energy correction in section 4, the Rayleigh- 
-Schroedinger perturbation theory is applied. Moreover the divergence types of the 
second order energy shifts are proved to differ from those computed by Brueckner 
(1958). 

Quite recently, the Green function method of Schwinger (1951) and of Gell- 
-Mann and Low (1951), has been successfully applied to problems of quantum mechanics 
Hubbard (1957b), (1958c), Galitskif and Migdal (1958), Klein and Prange (1958), 
Kanazawa and Watabe (1960), Kanazawa, Misawa and Fujita (1960). 

The problem of the Green function method in application to an electron gas of 
metallic density will be dealt with in a subsequent independent paper. 

(23) 
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2. The canonical transformation to independent collective modes 


The notation throughout the present part is the same as in part I hence no expla- 


nation will be given here. 
The Hamiltonian of the system, as derived in part I of this paper (I 3.8), is 


Hz YE (Ry MERE + 


|k|>k, 
RZ 
e TT 3 ^ N ik - r. 
oy V DĄ > qr &r (Pi — hkj2) ek "i — 
m 
=1 |k|< 
2ze? RER): r; 
NE PrP-» + © qk QE) + =P S SE dk qe ER" Ep EET)" i + 
k|<k > i=1 |k|<k, 
PTS 
k+Rk' 40 
N' » p, 
27e? ; N Ea ce 1 
W my 70 D > Lade ROUE ee V ko 
i=1 |k|Zk, LESS 
NASE | Sie G | dr dr” a _ 4a N'e? (2.1) 
= A. = —. CAE 
V por .2F9 fi ror : ml 


The frequency of plasma oscillations w, is, of course, not the same as in Bohm 
and Pines' paper (1953b), since here the number of individual electrons is smaller. 
Denoting the Langmuir frequency by c, we have 


4r Ne? N’ Ÿ 
2 = —- 2 ies BS 
SE mV N e. ( 5 P ) CH 
wherein 
Ar Ne? 
wL = mp (2.3) 


and Ÿ and f are given by eqs. (E 3.11) and (I 3.15), respectively. Thus, a certain 
shift due to excitation of electron-hole pairs is seen to occur in the plasma fre- 
quency wyr. 


Our main task now consists in transforming the Hamiltonian (2.1), so as to elimi- 
nate the interaction term 


^ e Ac CS. GUMMI ; 
Hines = 1/7 >, à qu Če (Bi — hlej2) etn. (24) 
i=1 |k|<k, 


The perturbation-theoretic treatment of this problem is due to Bohm and Pines 


2) 


(1953b). The generating function of the transformation is known to be of the form 


N' 


^ h , ar 
S = ELT. hy: k | Ek * (b; hk/2) ei: "iq, — 
m wr — k-pilm + hk?/2m 
i=1 |k|<k 
— qt e-iker, ër (Pi — hk/2) e 27 
x wr — k: film + hk?/2m |’ E Vhoz - k? ^ 2) 
We transform the Hamiltonian (2.1) as follows: 
Dan = eh H eih, (2.6) 
€, is the dispersion frequency satisfying the relation 
| Ame? NS 1 pe 
|. mV = (ox — k - Dim)? —h?k?/4m? EP 
or 
1 Ree 
Op = Op joe LC ET PT D D Pi + TT +. 3r (2.8) 


The right hand side of eq. (2.6) may be expanded in a H, S commutator series. Each 
n-fold commutator 


[H, Ś], = [... [[LH, 5], 125], BŚ] ..., 1^ S] (2.9) 


is proportional to the operator (k : p,/mo,)". The mean of the operator (k : P,/mo,)? 


1 2 1 
g? = < (k - p/mo)?> ~ <(k- fymo,)?> = 0.5426 E x (2.10) 
kapa y 
(+5) 
a = k,lkp (2.11) 


will be termed the coupling constant for the electron-plasmon interaction. 
Applying the random phase approximation (Bohm and Pines 1952a)), and the 
canonical transformation with the generating function (2.5), we obtain 


kom —r;l 


er l 2 |, sin T 
E EET 2 uni id 
Tew = > Im* + e , j (1 E | dt > | Els 
$=1 i<j 0 
^ il 2a N' e? il 2a N e? 
+ ho; 3 ak | 3) V 2: po DES my qo — 
[I|  & k|<k, 


al N'2e2 [dr s e? dr ar IDE gr (2.12) 
V 2 ya 
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wherein the effective mass m* is given by 

m V ke | in 3 
apa) (1-3 | coe 
A detailed derivation of eq. (2.12) is to be found in the paper by Bohm and Pines 


(1953b). 


Since we investigate the ground state of the system only, we can transform the 
oscillator part of the Hamiltonian (2.12). Using the plasmon part of the wave function 


(I 3.9), we have 
<s> = 6% > hoy (ak a4 + 1/2) D, = 4 D hwg 
| 


IKI<k, ki<k, 
or, with respect to eq. (2.8), 
ake 
^ ho k? D? Le h V p 
— moż Lei 2m m? w; 2 L(2zx)* . 
|k|« c 1 |k|<k, 
N’ A 
1 BV Rat 
| t a. 12.1 a. [4 
LAN: mo? > 2m (2x)° J ae 8m? w2 (22)? EE | 
= |k|<k, : LESS 


ho, | V An ,s 1 OF V 4a he V 4x 
EG | 3g he ty DN ke + : midi ia 
2 |8n? 3 N' moj 2m 8a? 5 8m? wh 823 7 ^ 


i=1 


ho, ya On 9 \ 
= anf de !vy3 ; 
5 N’ a? + 10 Ü | 5 | 


56 2 


a 


/30 p [929 h a? 1 e? 
| : 2s PE (2.14) 


We may now delete the fifth, sixth and seventh terms in eq. (2.12) as they only serve 
to cancel the Coulomb integral of the second term. By eq. (2.14), we have 


klr;—r;l 


NK „A 
D P? à ] 2 sin T 
"ue 2 2x pe 3 Iri — r; ( 7 | SĘ F 
i-1 imd J u 
0 
APSR Y (2.15) 
wherein 
IRE d V3 ono ^ a5 l 
x L- Far b o (n? Gw 
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V3 , [9x9 \" a? 1 
Si M Lr wer PRZ Ex) 
1 uc zo 2 


We have transferred one part of the long-range Coulomb energy (2.14) to the kinetic 
energy of the individual electrons. This step implies the dependence of the effective 
mass on the spacing parameter r, as in eq. (2.16). 


3. The correlation energy 


All successive computations are the same as in part I; hence, we restrict ourselves 
to given the final result for the non-magnetic state only: 


p — [2200 | 1-ie® , 0.3190-a8 0.9163 1 
Eoque marne] n Qe 
1.2217 a 0.8661 ^ a? 
AZ I da Fe 
APR ER A n (ie rh (cda tr 
0.2098 - a? 
1 7 TT | : (3.1) 


The first term in eq. (3.1) represents the mean kinetic energy per one particle, the 
second — the exchange energy of short-range interactions, the third — the self- 
-energy of long-range interactions, the fourth — the Coulomb long-range energy 
and the fifth — the correction due to dispersion of the oscillation frequency. Here, 
the designation *kinetic" energy is employed in a somewhat formal manner, since 
the first term in eq. (3.1) includes, in addition to the pure kinetic energy, certain 
interactions, such as i.a. a part of the Coulomb long-range energy. One part of the 
exchange energy (the second term in square brackets) cancels the long-range self- 
-energy (the third term in eq. (3.1)). 

We now proceed to minimalize eq. (3.1) with respect to æ. The corresponding 
values of & are computed numerically for r, = 1, 2, 3, 4, 5, 6, Table I. The remaining 
points for alkali metals are computed from the interpolation curve 


a = a + da'r, + ay: r? + ag: r + az" r; + ag". 62) 
Quite similarly the correlation energy 
w 9 
Eu EE, (3.3) 


| 2.2099 0.9163 


r2 r 
5 $ 


E Ry, i (3.4) 
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1,2, 3,4, 5, 6 directly (numerically), whereas, for the alkali 


is evalueted for r, — 
metals, it has been obtained from the interpolation curve 


E" 


COIT 


= b, + 0.0622 In r, 4- b, * r, + bg te + by: 1, + bz * fr + bg’ JRE (3.5) 


Putting by turns r, = 1, 2, 3, 4, 5, 6 in eq. (3.5) and denoting the corresponding 


values of æ and E by a(? and Ef, respectively, we get, after some amount of | 


algebra, 
(LIT Picks Ts |e E es | eee TOP a | 
| 172 12428 aha | by | |? E EC — 0.0431 
11444441 fo fey | E — 0.0862 
| Et E b, l| Fr — 0,1001 
176. 6400| l) [Eg —o114] 
(in, Sete UNS 
|. 10 6 24 6 5! 
117 461 137 19 1 
pi 4 24 24 24 4! 
127 121 3 1 
3 SP = : Ex im repe Lic 
jp e can 3 d 12 2 12 em 
di : 307 107 1 i 
rs er 12 2 12 
A 27 65 95 2 1 
bac 6 24 CE ae 
^ 137 15 No EL E 
Fe 60 8 24 8 51 


wherein A is the determinant of the square matrix in eq. (3.6) and A — its minor 
(i and k stand for the lines and columns to be deleted in the determinant A), 
Numerical computation yields: 


Table I 


zę l 
a 0.4887 
E, Ry 1.2339 
E Ry 1.2936 
EM. Ry| —0.0597 


2 3 4 5 6 
0.5748 0.6454 0.7034 0.7527 0.7975 
0.0622 —0.0856 — 0.1133 — 0.1151 — 0.1097 
0.0943 — 0.0599 — 0.0910 — 0.0949 — 0.0913 

— 0.0321 — 0.0257 — 0.0223 — 0.0202 — 0.0184 


29 


Since 
6 : 
„ At 
que a e (3.8) 
| 
^ 54 
4 1] 
b, = D ej p (3.9) 


we may now make use of the values in Table I, obtaining 
« = 0.3859 + $ 0.66591 : r, + 4; (—0.1891 - r; — 0.0103 - r? + 0.0031 - r$) — 

— 1% 0.0007 : 7°, lm <0, (3.10) 
Ecórr — 0.0747 + 0.0622 In r, + 3 (0.13889 - r, — 0.1265 : r? + 0.0374 +73 — 

— 0.00505 : r$ + 0.00026 : r5) Ry, DOT 6: Gall) 


From (3.11) and (1 3.54) the constant in the correlation energy eq. (3.11), is seen to 
be by approximately 20% larger than that computed by Gell-Mann and Brueckner 
(1957). 

Table II lists the interpolated values for the correlation energy eq. (3.11) for 
some alkali metals. 


Table II 
Metal Li Na K Rb Cs 
r, 3,22 3,96 4,87 5,18 5,57 
0,6591 0,7013 0,7466 0,7610 0,7785 
E, Ry — 0,0956 —0,1128 —0,1153 —0,1142 —0,1124 
E Ry — 0,0715 —0,0905 —0,0950 —0,0945 —0,0933 
EM. Ry —0,0241 —0,0223 —0,0203 —0,0197 —0,0191 


The values of the correlation energy, as evaluated either with the help of eq. (3.10) 
or by eq. (3.11), are consistent within 3%. 


4. Second order correction to the correlation energy 


Let us now investigate the magnitude of the second order energy shift. For this 
purpose we utilize the Rayleigh-Schroedinger perturbation theory (Brueckner 
(1958)) and the wave equation 


EK E (4.1) 
NG. sa N’ 
^ b; hh? ke 
og 2; p Cya FU Dt 2 
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with the wave function V given by (13.9). According to eq. (2.15) the Coulomb short - 


-range interaction energy is of the form ' 
klr;-r;l 
e? 2 sin T 4r e? ekir 43 
h = hy, huj © dt = — RB , (4.3) 
5 Ir; — r;| c T I 
i<j 0 |k|>k, 


whence we obtain the following energy correction: 


1 1 
EN EA I de he 4.4) 
| b ) A ROT 
with 
1 E l — Po P, É © when m V. (4.5) 
b. E. fo 0 otherwise. 


b, is the projection operator on the unperturbed state V^. By (4.3) and (4.4), we have 
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In Gell-Mann and Brueckner’s paper (1957) all terms in the second and third sum 
vanish, since both authors put 


N 
D I qu) galt) pet (4.7) 
i=1 

It is manifest that the wave function in the form of eq. (4.7) does not obey the Pauli 
principle. Thus, in their theory, the correlation effects due to the Pauli principle are 
neglected. In the present author’s computation, the aforementioned terms are non- 
-zero; albeit the evaluation of these terms is excessively tedious, and they will be 
omitted in the present paper. This yields 


4x e? | 2 1 ai LE 
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Insertion of eq. (I 3.9) into eq. (4.9a) gives 
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With respect to eq. (4.9b), the expression (4.8) becomes 
£y = eW + ap, (4.10a) 
where 
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Averageing the two second order energy shifts (4.10b) and (4.10c) over all values of 
pe wave vectors k, and k, within the Fermi sphere (I 3.12), we obtain (per one indi- 
vidual particle), 


0.6079 i 
= 0000 — ——— — 0.5120 : x — 0.1013 - aż + 0.1230 -a° + E 0.0853 + 
0.3843 0.6321 0.3226 + a 
T — | +24 — | 0.0161 + — a] <a = 
BUF: ca i | prex d TESTER 
0.1281 - a 
€ ~ 24 170 cw! 
m [020% a — 0.0034 - a3 + 0.0338 - a z SIM nal Ry (4.1la) 
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0.0169 In? a | 1.0845 — 0.1695 - a + zc e] In (2 — a) 


— 0.0169 In? (2 — a) Ry. (4.11b) 


Both quantities are evaluated in the Appendix. 
Adding e@ and e$”, we have 
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It can be seen from eq. (4.11a) that &&? has a divergence (x — 0) of type a1, whereas, 
according to Brueckner (1958), it is of logarithmical type one, 
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This discrepaney seems to be due to Brueckner's approximate integration as opposed 
b 

to the computation given in the Appendix to the present paper. As to the term eO, 

eq. (4.11b), it diverges as In?@, whereas Brueckner found it to be regular at œ = 0. 


The second order energy correction eq. (4.12) depends upon the spacing parame- 
kk 


eq. (2.16). 


ter r, due to the presence of the quotient 


In order to compute & we utilize the values for a listed in Table I. This yields 
the following Table III and the Fig. 1. 


Table III 
r 1 2 3 4 5 6 
ela) Ry — 0.5699 — 0.2920 — 0.2583 — 0.2337 
e@) Ry — 0.0109 —0.0107 —0.0114 —0.0119 
e, Ry — 0.5808 — 0.3027 — 0.2697 — 0.2456 
5 6 


0.10 


% 


025 
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Fig. 1 


It is now clear that the second order energy correction (4.12), if added to the 
correlation energy, lowers the latter many times. This is for the most part due to 
e$), It may be hoped that the successive higher order shifts would reduce the enormous 
contribution of ef, but it is not possible to verify this conclusion, since the Rayleigh- 
Schroedinger perturbation series for r, 1 is divergent. Therefore we choose the fol- 
lowing procedure: in part I we made use of thefact that the longitudinal vector potential 
is defined to a constant and we introduced its zero component go, which we put 

k.|r=r'| 


2n N'e? N (N' — z dr dr’ sin 7 
mV gij aV? ls — ud f gr ar (Som) 


so as to obtain zero electrostatic energy. Since the term &® is of Coulomb integral 
type structure, we may now extend the definition of go 


klr=r'| 


2x N'e* , N (Ni 1) dr dr’ sin T 
my 40 ~ aV? IE Zr] Í dt = muet, (4.13) 
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z E B b . . . . 
Thus, only the second correction e? remains, which we incorporate into the correla- 
tion energy. This yields 


EX, = {—0.0921 + 0.0622 - In r, + $ 0.34373 +r, + & (—1.1449 + r? + 
+ 0.3328 + r3 — 0.0445 - 7) + 3 0.00114 5) Ry. (4.14) 


The constant in the correlation energy is now by nearly 4%, larger than that of Gell- 


-Mann and Brueckner (1957). 


5. Conclusions 


In the first part of the present paper, a collective approach without subsidiary 
conditions is proposed. It was assumed that the collective oscillations of the charge 
density are excited by the colliding electrons and the electron-hole pairs appearing 
thereafter. Both phenomena — the electron-hole pairs and the collective oscilla- 
tions related to them — were described with the help of the effective longitudinal 
field. The theory thus formulated allowed to calculate the magnitude of the correlation 
energy to the first order. 

In the present paper, we effected the isolation of the plasmons and computed 
the second order correction to the correlation energy. Up to the first order, the constant 
in the correlation energy is approximately 20%, larger than the same quantity obtained 
by Gell-Mann and Brueckner. Together with part of the second order shift, the constant 
becomes larger by about 4% than that calculated by these authors. 

As we have already mentioned, Gell-Mann and Brueckner applied approximate 
integration when averageing over the Fermi distribution, considering the probability 
for exciting the electrons from the neighbourhood of the Fermi surface as being 1 and 
the probability of raising electrons from the depth of the Fermi sphere onto external 
orbitals as being zero. Such approximation gives rise to some error in the correlation 
energy constant. 

The method of determining the correlation energy constant (interpolation curve 
for the correlation energy) used throughout the present paper involves a degree of 
uncertainty residing in the choice of the expansion in r,. No such complete expansion 
is available, neither is one given in the Gell-Mann and Brueckner theory, as the latter 
retains linked cluster subseries only, whereas the processes called “non-ring”, includ- 
ing correlation effects due to the Pauli Principle, are neglected. 

The author wishes to thank Professor Dr S. Szczeniowski for his interest through- 


out the present investigation. 


Appendix 


We first prove the correctness of eq. (4.11a). For this purpose we write ein the 
form of eq. (4.10b). 
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Let us, for the time being, write 


k 


= = 0 
p 70-302, 


* where 2 is the unit vector in the direction of the z-axis. We have 
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We now average the expression (A.1) over the Fermi distribution (13.12) for à = 4, 
obtaining 
V 1 z 2 1 ] ke +k, 
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eq. (A.2) transforms to 


ke + 0 b 
x ln dR. A.A 
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Fig. 2 Allowed region We evaluate the second integral of eq. (A.4) in cylindrical 
(shaded) for R integration coordinates, 
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Reverting to eq. (A.4), we have 
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All integrations in eq. (A.6) are elementary except one integral, namely 
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The last two integrals are readily computed by expanding In 
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(Jahnke and Emde 1959). 


In order to complete our computation, we still have to expand In E za | and 
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In [1 — € | in eq. (A.6) (we retain only the terms up to a5) and to insert the expression 


obtained into eq. (4.10b). The result is eq. (4.113). 
We proceed now to the evaluation of the energy correction (4.11b). We have 
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We assume the k-vector to be directed along the z-axis. There are three stages of 


integration illustrated below 


Fig. 3 Allowed regions (shaded) for a integration for: a) k,<k<2Żkp— ko 
b) 2kg — k, <k<2kp+k,, c) 2kp+ k, <k < + oo. 
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Since the contribution from small K are the most essential for the integral (4.10€), 
we restrict our computation to the case a). The integration now consists in taking 


the integral over the sphere of radius 2%, minus the integration within the k,-radius 
sphere. We have then 
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In the term (A.7e), we substitute K -- a = b obtaining 
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Since b < k, < k, we may expand the subintegral expressions in eq. (A.7g) in power 
series in b/k. Thus, we have 
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With respect to eqs. (A.7f) and (A.7h), we obtain 
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we obtain, per one 


Inserting eq. (A.7i) into eq. (4.10c) and multiplying by £ 


particle, eq. (4.11b). 
It may be readily verified that a contribution from the case c) of Fig. 3 is strongly 


convergent and yields a small correction only. 
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The probability function of occupying the states with definite energy in the conduction 
band after an ionization is calculated. The valence band is assumed of finite width. Several 
cases of the primary energy and of the effective mass of the valence band are examined. 


The electron impact ionization plays an important role in such processes occuring 
in insulators and semiconductors as: (1) the electric breakdown in insulators (see 
Franz 1952), (2) the interaction between radiation and insulators or semiconductors 
(see Antoncik 1957), (3) the processes occuring in semiconductors under the influence 
of electron bombardment (see Zaręba 1959). 

The electron impact ionization is a process in which a sufficiently energetic 
electron from the conduction band, or alternatively an electron coming from the 
outside of the crystal, suffers a “collision” with an electron from the valence band, 
so that the latter goes over to the conduction band while the primary electron remains 
also in the conduction band. 

The theoretical treatment of the problem of the electron impact ionization was 
given in Tewordt's paper (1954) and later it was summarized and generalized by 
Franz (1956). 

We shall follow here the definitions and the notation used in Franz's paper. 

Franz has given the expressions for the probabilities of the electron transitions 
in the process of electron impact ionization. Their derivation was based on some 
simplifying assumptions. One of the fundamental assumption is that the electron 
impact ionization is considered only in that region of the wave vector space which 
is close to the ionization threshold. 

Franz has also given an expression for the probability of occupying the states 
with definite energy in the conduction band after the process of ionization. This 
expression was given only for the case when the valence band is treated as a discrete 
level with a width equal to zero. It follows from this expression that the most prob- 


. (43) 
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able distribution of the surplus energy is the equal distribution betweea the two 
“colliding” electrons. 

In this paper the probability of occupying the states with definite energy in the 
conduction band is calculated for the case of a non-zero width of the valence band. 

We shall perform our calculations for the case when the ionization is caused 
by the electrons in the conduction band with a small energy E(k,), namely J < E(k;) < 
< 2J, where k, is the wave vector of the primary electron. The zero of the E(k,) is 
taken at the bottom of the conduction band. Franz's assumption is followed that the 
electrons in the conduction band behave as free. J is the energy gap between the 
valence and the conduction bands. For the energy E(k,) within the specified limits 
the ionization takes only place in the region of energy close to the threshold of ioniza- 
tion and that takes place near the top of the valence band. In this case the assumption 
of a parabolic form of the valence band energy is a reasonable assumption. Furthermore 
in this case all the approximations made by Franz are well satisfied. 

Therefore, the following expression for the energy of the valence band is assumed 


h? 


2m* 


E(k) = —J— 5 M. (1) 
where k, is the wave vector and m* is the effective mass of an electron in the valence 
band. 

The following situation is considered: an electron from the conduction band in 
the state determined by the wave vector k, interacts with an electron from the valence 
band in the state ky. After the “collision” there are two electrons in the conduction 
band in the states determined by Kj, and ką. In the process the laws of conservation 
of energy and momentum (in our case the latter is equivalent to a relation between 
the wave vectors) must be satisfied: 


ki + ką = k, + k, + k, = 2ky, (2a) 
E(k,) + E(k) = E(k,) + E(k,). (2b) 


Here k, are the vectors of the reciprocal lattice. These two laws are fulfilled simulta- 
neously if 


l 2m d 
2 : 
R= >) ei de E (ky) | — ki 70, (3) 
where Ko is a function of (k, + k,). From the condition K?> 0, we can determine 
in the valence band the region of the vectors (ką + K,) (or alternatively ko, which 
we shall use later) from which, for a fixed k}, ionization is possible. 

The probability of occupying the state with energy E(Kj) when the primary 
electron has the wave vector k, is given by the expression 


2x 
| V 1 
P(k,, E(ki)) = ee Ji dz ko | dp zg P(k,, ko O), (4) 
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where P(k,,k,, 2) = DD (o Hor E k,) has been replaced by (2k, — k,) and 
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P, (ką, ką, Q) F: Vh3 a(k, + k;)| K lei + 
Be = : a) 
BRIE (E, — E, — Ej (ke, — k, TK e) 


a(k, + k,) is the Fourier component of the wave function of the valence electron. 
With the approximations employed by Franz, the final expression, which will 
be used for the calculation of the probability function, has the form 


P(k,, E(ki)) = C fl 3 (6) 
0 


RAM, 
4m? ef (UE Qa ||| TH F 
z?h? (1 + u)* ki , (7) 


where 


2 


and u = m/m*. 

In order to calculate the above integrals, we must know the region of integration 
in the Ky space. The region of Ky vectors allowed by the conservation laws is obtained 
from the condition K?> 0. 

The primary electron with the fixed wave vector k, transfers its energy to an 
arbitrary electron with the wave vtecor k, (Ky substitutes k, + k,) from the allowed 
region. However, it is easy to verify that an ionized electron with a definite energy 
E (Kj) cannot be ejected from every state in the allowed region of k, but only from certain 
states. In other words, the ionized electron can attain the energy E(Kj) by a “colli- 
sion” of the electron in the state k, with any electron from a fraction of the allowed 
region of k, only. Thus we shall have different regions of integration for different E (Kj). 

Turning our attention to the function P(k,, E(k;))}, we must discuss its 
dependence on one variable E(K;) and two parameters. The parameters are: (1) 
The momentum or the wave vector k, of the primary electron. Due to the spherical 
symmetry of the surfaces of constant energy for both the conduction and valence 
bands the dependence of P(k,, E(k;)) only on the magnitude of k, or on the energy 
E(k;) can be considered; for every direction in k-space the situation is the same. 
(2) The parameter u or m* — the effective mass near the top of the valence band. 

We can go over now to the discussion of the region of integration. Firstly from 
the condition K2>0 we shall find the whole region of the allowed k, and in fact the 
region of 2k,, because this latter vector appears in the relation (2a). 

'The solution of the equation K? — 0 gives the boundary surface of this region. 


Solving this equation, we obtain 


7:91 - 
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fig : : . 
= ką is the smallest. energy of the conduction electron. for which the impact 


2m 


ionization begins. It depends on the size of the energy gap and on the effective mass: 


a = 2 
HR Er 


D Ml PE 9) 


We see now that the region of the allowed Żk, is bounded by a sphere, the centre 
of which is shifted through the vector Żuk,/(1 + 2u) relative to the origin of the k; 
vector. (This origin is taken as the origin of the co-ordinate system in k-space). The 
location of this sphere relative to the origin of the k, vector is important. The origin 
can be either inside or outside of the sphere. It is also important whether the end 
of k, is inside or outside of the sphere. The location inside the sphere indicates 
that ką = 0 is possible. This means that the ionization can take place from the highest 
level of the valence band. It is easily seen from (8) that for the energy of the primary 
electron considered here, namely J < E(k,) < 2J the ionization from the highest 
level is impossible. In our case the end of the k;, vector is always situated outside of 
the sphere. 

The following cases may occur: 


I. Shown in Fig. 1. 


| wide == (2 Kymin “| 
| P, 


It takes place for the following m and K,: 
O<u<ż H«H«JI(-—y) 


ENIE <<). 


— m 7 
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For the magnitude of 2k) we have the following analytical expressions 


2ko.. (6) = x ; E cos f AM E [AŻ (1 — p) E 


Es = 3 ky {cos B + cos? B — cos? Pr}. (10) 


The angle f varies here from 0 to f, and 2ky, refer to the two parts of the sphere in- 


dicated in Fig. 1. 


34 V - 0-09 I. (11) 


The minimum value of |2 K,| on the sphere is attained for ß = Oand equals2k,_ (8 = 0) 


d unosi ZU Ro FIDA 9 
Cian = Den — V1 — cos? fr), (12a) 


and the maximum value is 2k,, (B = 0) 


Pho)max = zy = 37 ky (1 + V1 — cos? fy). (12b) 
2ky is the value of |2k,| at B = fj. 
2 
PIE. ei łą cos Br. (13) 


II. Shown in Fig. 2. 


| 
| 
| 


| 
= (2 ky max 


Fig. 2 


Case II occurs for u < $ and JI ED ZE 2. 
For 2k, we have the following analytical expressions 


Hoo (P) = pzy [t ]/ n ki cos? B + mud [ki (1 — u) — 1 
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=; = + “co B+ + V/cos? B + By? yt. (14) 
|e 2u DN 

- ABN " i 15 

24 ki POM E 


Here the angle f varies from 0 to z. 
The minimum value of |2k,| on the surface of the sphere is attained at f = 


2 Piece dora 
konia = z MH VIE) (16a) 
and the maximum value is at f = 0 
2 ue 
Gau = zy h (L + VI y. (16b) 
l + 2u 


We now look for the region circumscribed by the ends of the vector (k; + k3) 
taking into account only the law of conservation of energy. The primary electron can 
ionize the electrons from different allowed energy levels in the valence band. Hence 
the total possible energy of this pair of electrons before the ionization process 


6 =E(k,) + E(k,) 
changes its value according to the energy level from which the electron is ionized. 


o o 
<6 Si aes 


o 
O min 
where 
e 


== = E(k) T = ky) = E (ky) b E(|Ka|min) 


Pg onse (17a) 
Ó min E (ky) AXE "min (kə) + E (Rh) 
2 
zr EM. (17b) 


and (Żko)max: (2Ko)min are given by (12a,b) and (16a,b). In the further calculations 
zur P 
pe ©: 

In order to fulfill the law of conservation of energy in the ionization process the 
vectors ką and ką must be such that 


we shall put é = 


Gmin 5 > E (ky) Y E (ka) < ? max (18) 


Since the ionization is possible from different energy levels in the valence band, 
we can have for a fixed E(k;) such E(k;) that the relation (18) is satisfied. 


E(k;) can vary from 0 to 6 max: At E(k,) = 0, 
dwa PH ' (ką) > a Ó min? at E (ky) = ZP, E (kz) = 0. 
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The region circumscribed by the vectors (Kj + kj) can be easily found from 
the relation (18). The following situations take place: 


1.0 <E(kı) <4 Onin. 


The ends of (ki + ks) give an excavated sphere with 


i See rare =e /2 
internal radius r,;, = V @min — ki 


y 
external radius R nax = IA = kę + ką ; 
2 i DR E (ki) < +6 max * 


The ends of (k, + ks) give a sphere with the radius 


R ax uc ven wu Ky 3x łą 9 
| ARNE ETA 


The ends of (k; + ks) give an excavated sphere with 


Ds = = 


external radius R, = ką oe N — kr 


internal radius R 


Since the vectors (ki + ks) must simultaneously fulfill the law of conservation 
of momentum ki + k, = 2k,, the region of integration for fixed E (ki) is given 
as a common part of the spheres discussed above and the sphere of allowed 2ko. 

The integral in expression (6) can now be evaluated. A more detailed discussion 
reveals that the regions of integration depend on the values of the ratio © min/Ô max: 


This ratio is a function of u and k. 
We consider first the regions of integration for the lst case above. There are 


four different limits © ,:,/@ max in which the integration regions are different. 


(1.2) 


| 


bole 


) = Ómin/Ć maz < 
Í < min/max < 3 (3 — V5) | (1.3) 
OE © min! max < (1.4) 


bole 
— 
W 
| 
= 
au 


Bie 


Now we enumerate for what values of wand k, the © nin’ max lie in the specified ranges. 
These ranges are determined by the characteristic parameters, that is the values of 


k, for which 
On| © thas T 2 > (3 y5) > : 
D min! max 5 3 for 


1+ y — 2u* (/2 — ne” 


Il D = 
Ita 2j (/8 — 2 
P dni mar PE i for . 
Lane 1+ 2u Tr 


1 PE Lao de 


Notice that x? < x? < x3. 
We can now write down the values of u and k; for different ratios 6434/6 max and 
therefore different regions of integration. 


0SK<; 
kż < < 22 AE RK mu 
xd? 4 (3—V5) < C i| max < à 
00 CR Gd 1 < min/max <3 (3 — V5) 
x < k <J/(1 — à) 0 < Omin/O max <} 
302 23 
B <k < zę 1 < Śmia/Ómsx € 1 
x < 12 < 8 1 (3—V/5) < Omin/O mex € À 
x3 < ki <a t < Omin/6 mex «1 (3 — V5) 
x3 <k? <2J 0 — Omin/Omax < 4 
s<Ku<1l+3/5 2.118 
k? < kè o EO N EA 
x < ki i 1 (3—y5) < @min/Omex € 
id < k? < 2] 4 < C nin[Ó max <4 (3 — V5) 
1 +4V5< u<} + V2 ~ 2.91 | 
H«HexX = Onin Sines 1 
x <k < 2] 4(3—V5) < Smin/Omax € à 


In the second case there are two different limits 6 min/max in which the regions 
of integration are different: 


1<6 = (11.1) 
ORCL lee EE (11.2) 


Here Cain © mex = à for 


TELS A 
kj = xh = c J. 


1+ p — 2u2 (1+ V2) 


We can now write down the values of u and k, for the different ratios & minl 
in max * 


0« pu <$— y2 ~ 0.0858 


JIQ — p) < E «x 0 < Ó min/Ó max < 3 

E wid 5 Omin/Ó max < 1 
Rope 

Tl uc e 2] 0 < ó min|Ó max < 3. 


The determination of the integration regions is shown in detail in Fig. 3 for the 
particular case (1.1). There we exhibit the dependence of R, R and Tin on 


12 Ko, 
Wa Scenes Trece te TEE (mik 
-—— k 
E maz va 
€ min 
| 
| 
| 
| 
| 
! 
| 
| 
qeu pcm 
VEminl2 VE RZA (ZMR 
Fig. 3 
k; and the situation of the curves oec SRL „ (ką) and ry „ (ki) relative to the 


sphere of the allowed 2k, The figure exhibits ab the Side of the vector 


2k, without taking into account the angle p. 


We can find from the figure the values of k, for which Ry, (ky) < kou 


that is the spheres considered have no common region. Let us denote by ką ,_ and 
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ką 14 the points of intersection of the curve Ry, (Kj) with the line (2k9)min« We 
notice at once that for the values of ky in the iterato 


, ios 


, , , [2] 
0 < ky < ką E — and ką a < ką = | O max 
. . B . u , 
the region of the k, integration is zero. Therefore for these values of A, the pro- 


bability function P(k,, E(k,)) vanishes. 


+ 1 e 
The values of k,,,, are given by 


, l = = = 
las = y2 Ka Æ ys. zę & min ) 


In Fig. 4 the k, integration region is shown for the case (2ko)min < Rmax (k1) < 2k. 
Considering Fig. 3 and denoting there by ka. , Aj. the intersection points of the 


curve Rmax (1) with the line 2kg, we see that the integration region of the shape 


max 


Fig. 4 


shown in Fig. 4 appears if the ką values lie in the intervals E rek ki < k12_ and 
Eu Shy A 1+- The values of ki a4 are given by 


bia. = y2 (Va, +VE max — 61), where 6} = pee “Can 


[he integration limits resulting from Fig. 4 are 


eB ep WER ER 


max * 


ne Sz fy depends on Kj. The cos s Bo i is obtained from the equation 2ko- (fq) = 
zm iis (ki) as 


—— l 4 24V ,, 
cos Bo TER Il 2; Rhax + ki cos? B |. (19) 
In Fig. 3 we have still to consider the interval of the ką values ky, vidéo MEA m 
In this interval Rx (Kj) > 2ko and it attains its maximal value Roz = 5 (2ko) max 
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at ky = /6,,,/2. The region of the k integration for this interval of £j is shown in 
Pie 5 
We have here the following limits of integration for f and 2k, 
0 «p < Pi 2k5- < 2ko i uus 


= << 
fı € p <8, 2kę_ « 2k < 2k 


ae , 


The angle f, is obtained from the equation 25, , (61) = Rmax (Kı) and is given 


by expression identical with (19). e 


In an analogical way, using figures similar to Fig. 3,4,5, we can determine expli- 
citly the regions of the k, integration for the other cases. 

Performing the integration for all these cases the expressions for the funetion 
P(k,, E(k;)) are obtained. 

Let us now collect all the characteristic values appearing in the determination 
of the regions of the ky integration for all the cases: 


J — — 
ki1+ p V2 Vas + pet. ae Ó min); 


, 
kia. 


Ve, + jr i 1), 


l 
SDE y2 ( 
kis = V ómar? + Gr. 


For the I-st case 


2 
E C CS: 1 + 24 2 1210: | 
cos By = cos p, = zw EUN E Il 2u Rinax + hp cos? fi]. 


cos fly is obtained from the relation Rmax = 2ko_(Po) and 
cos f, from the relation Rmax = 2k, (f,). 


2 
© Eum. nus | 
cos fj, = cos ps = EXEAT. Il zu in + RL COS" Bj 
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cos fj is obtained from the relation Ry, = 2k_(8ą) and 


cos ß, from the relation Rmin 


= 2ko+ (83). 


In the case II we have in place of cos? fj the quantity (—y?) in the expressions 


for cos fg and cos fo. 


Let us now write down the function P(k,. E(K;)) for all the cases mentioned 


above. 


Intervals of ki 


P(k,, E(k) 


, Fi 
lio , ke) 
, 


, 
ki 24 > ky 44) 


( 
( 
( 
Pen ? "EN 
( 


(kyo , ki o) 
(sx , k 14.) 
Can : IE 
(0, kia) 

(E a , ką i4) 


, , 
Uu > hy 24) 


(Has V dan) 
(0, b +) 
(CPR , ką s) 


, , 
s , ki » 


LE V Emax) 
NE 


mesic EIS 


AU 
1,1+> V © max 


V Sial 2, VE mul?) 
Von 


0 (1.1) 


CA; (PE ME cos Po) 
CLA, pede cos Pi) + As 
C. A, (Rak COS Do) 


(cos f4)] 


CA, (Rs cos Bo) (1.2) 
CTA, (Ry. cos By) + 
C. A, (Koax, COS Bo) 


43 (cos Pı)] 


C[As (Rmax Rmin» COS Pa) + Ag (Raz COS Bo, COS Pa)] 

CA (RS T =A) (1.3) 
CIA (Renae cos By) + As (cos Bj) 

C[A5 im " + As (Rao Pig: £08 Ba) + 


+A ; COS By, cos Bal 
CL 1; i c min» COS pa) zh 1, ER cos Bo COS £2)] 
CTA, (Roux > COS By) + As (cos P3] (1.4) 


C[ A, (cos fi) E Az Bier Ii COS P>) d: 
+ tg (iy COS Bu COS Doll 
CTA (cos ba) + As f ie c NEWEST 
+ 45 n , COS Bj, cos f4)] 
[Az Raar Rain» COS By) + As (Ram. COS By, cos f3)] 
CIB; f sass "min > COS Bo) bi D, (CHARS cos Po)] (1.1) 
CIB; (Rmaxs cos Po) + B, (cos B,)] 
C|B, ME TR COS Bo) nx B; (Rs cos Ból 


C[B, (Rinax> Rmin> COS Bo) + Bg (Rain COS Po cos Bo)] 
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(0, ka) C.By (Ras ri) (11.9) 
EU) = GE tse Pains cos Po) + By (ass cos Pi] 

(Vus. Vsus)  CIB, Russ cos By) + By (cos Po] 

IE T UC RR ch) (Rag, tos P] 

CEPR [m C.Bg (Bax: Rein) 


The functions A, and B; are given in the following forms: 


A (Rmax cos Bo) vt Ę Ria (1 — cos Po) om wów ER ) łą E (1 = Ove RE 


— cos? f; (1 — cos fy) — 8 V (1 — cos? b)? + a V (cos? By — cos? m ||. 


4 z ——— 
Ay (cos fs) = 5 m l £x] H V (cost, — cost f 
A, diee Jimi: COS Ba) = AT „Alike oe TE (1 — COS p) 
A, (R cos Di. COS Bo) = z Em (cos B4 — cos f4) — "un ded Mi 
4 max» COS Pz, oj 4 max 2 1 Tee 


ES E (cos? B, — cos? B4) — cos? fj; (cos B4 — cos fy) — 


2 2518 2 —— Ó 
— 3 V(cos* By — cost)? + = V(cos? f, — cos? fi) | 


2 
2. 
À; (Fin; cos pi. COS Po) = ls d ki E (cos? By — cos? pa) = 


— cos? fj; (cos B4 — cos fs) + > V (cos? B, — cos? Bj — 


— = V (cos? Pa — cos? po | = 2 Rin (cos f4 — cos z 
Be (Bmax; Amin, COS Bo) = (Ria — Thin) (1 — cos Bo) 


3 
B, (min; eos Bo) = % — an A (cos? Bo + 1) + y? (cos By + 1) + 


1 
+3 Yes AU tr — Y zyj — Thin (cos Bo + N) 
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l 
Bs (Rmax, COS Bo) el 4 Riaz (1 — cos Bo) 


2 
| NE 
B, (cos Bo) = x rza) ke B (cos? By + 1) + y? (cos By + 1) + 


2 2 VAL 5 
e + V(cos? By + y2)* — 3 ya + zy 


PEREZ 
B, (Rain; cos By, cos f3) Sr | im E T cos? By — cos? Po) + 


T 2 xe 
+ Y? (cos By — cos fa) + 3 V (cos? By + y3? — 


2 = 1 mo 
s V(cos? B, + zy =. Rhin (cos By — cos A 


l 
Bg (xs hon NS Ub (Roux w Ado). 


H / , . . . . 
The character of the function P(k,, E(k;)) was investigated in its dependence 
on E(k;) for several parameters u and kı which correspond to the different cases, 


of integration. 
We get only three characteristic, qualitatively different, shapes of the function 


P(k,, E(k;)), namely: 
(1) The function P(k,, E(k})) is for some values of E(k;) equal to zero. This case 


comes about for such u and ki < Ja — p< 2J, for which $ LO minal Oar EB 


100 


Ena f 2 ue / Ek A, 


Fig. 6 


1 1. . . ^f . . 
Fk) is confined to the interval (0, may): The function P(k,, E(k;)) is symmetrical 
with respect to its maximum. It reaches the maximum for Elki) = $,,/2. An 


example of such a function is shown in Fig. 6. It has been computed for u = 2, 


and ky = 1.297 |/J. 
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(2) The function P(k,, E(k})) is equal to 0 only for E(k,) = 0, & max- It is a symmetrical 
function with a maximum at E(k,) = =. This type of function is obtained 
for such u and ką, for which @min/Pmax = 3. An example computed for jmd 


and k, = 1.3452 VJ i is shown in Fig. 7. 
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Fig. 7 


(3) The function P(k,, E(k})) has a definite value, different from zero for E(K;) = 0. 
This case can be realized either for u and ki < Ja — u) < 2J such that 


OminiOmax € + or for p <4 and OJ ki JK 1— u). The maximum of 
P(k,, E(k})) is at E(Kj) =6,,,/2. The function P(k,, E(k;)) depends only 
slightly on Æ(k;) in the neighbourhood of the maximum. This means that the 


bs 


708x p(, 


4 Era 2 Enar E(k, ) 


Fig. 8 


probability is almost independent of E(k;) about 6,,/2 though it depends on 
E(k;) strongly about 0 and 6%. 
Two examples computed for u =0.1 are shown in Fig. 8 for kı = 1.0449 VI 


and in Fig. 9 for ky = 1.0452 y J. 
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Emaz/2 e ANGE, 


Fig. 9 


The expression for P(k,, E(K;)) has been given in Franz's paper (1956) only 
in the case of a vanishing width of the valence band. In order to get Franz's result 
from our computations, we must take a properly defined interval of E(k,) and consider 
the limiting case u > 0. 

We have to consider the following limits for u and ki 

0«uci-—y2 

x < ki <2J 
When u —> 0, x} > J and ką varies from J to 2J. In order to give the expression 
for P(k,, E(k;)) for u = 0 we must consider the case (II.2). 

When u — 0, then d, = nia = ó = Elk) — J, E = kii = 6/2 and for 

all the values Æ(k;) from 0 to 6 

P(Ie Elk) = Cum ODE RE 
R Am? et ak)? E 
EIL ką 2 


This expression agrees with that given by Franz. 


we have 


R 


mial) F 


TO. ps Ras) = 
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Academy of Science in Prague. The author would like to thank Dr E. Antonëik for 
suggesting the problem and helpful discussions. 


REFERENCES 
Antonéik, E. Czech. J. Phys., 7, 674 (1951). 
Franz, W., Z. Phys., 132, 285 (1952). 
Franz, W., Handb. d. Physik, Bd. XVII, Berlin 1956. 
Tewordt, L., Z. Phys., 138, 499 (1954). 
Zaręba, A., Bull. Acad. Polon. Sci. Cl. III, 7, 255 (1959). 


Vol XX (1961) ACTA PHYSICA POLONICA Pascal 


NOTE ON THE EFFECTS DUE TO PULSE-DELAY PHENOMENA 
INZCO, oe GM COUNTERS 


By W. Moscicx1 
Gdansk Laboratory, Institute of Nuclear Research, Polish Academy of Sciences 
(Received June 14, 1960) 


The possibility of using CO, as principal agent for filling G.-M. counters was 
proved by Brown and Miller (1947). Of many mixtures tested, CO, containing an 
admixture of several %% of CS, yielded the most satisfactory results when used for 
filling counters. The latter then presented a plateau length of several hundred volts 
and a slope of the order of 0.5—2% throughout the range of 50—100 V above 
the threshold Mościcki (1958). When measuring the Typ of 14C, Hawkings et 
al. (1949) determined the yield of counters of the type considered, which they found 
to amount to 96.2% in a counter of a diameter of 1.167 cm and to 98% in one of 
a diameter of 6.998 cm. 

However, it should be stated that Tjj, measurements carried out with these 
counters showed discrepancies, as compared with those of authors using other 
techniques, especially argon-alcohol filled counters. 

Direct comparison of the number of counts from the same sample as 
obtained with CO, + CS, counters and with argon-alcohol counters containing 
a small admixture of active CO», proved the former to present a yield of 80—85% 
of that of the latter type of counter. The reasons for this reduced yield remained 
unknown (Miller et al. 1950). 

When investigating the use of CO, + CS, counters in the MC dating method, 
the author noticed that the number of pulses of the background depended 
evidently on the time of blocking of the anticoincidence system by pulses of the 
shielding. Observation proved the number of additional background pulses AB to be 
given approximately by the distribution 


AB — 0.65 L e- 9.72 - 10^ (min) (1) 


with Z denoting the total number of pulses appearing in the counter, and 7 
(in minutes) the width of pulse of the a.c. channel. With respect to the low accuracy 
in the evaluation of the a.c. multivibrator gates width occurring in the investigation 


(59) 
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and the scarce statistics available, the accuracy of the constants in eq. (1) was 
estimated to be about 20%. 

The increase in number of background pulses when the a.c. pulse width under- 
went shortening was interpreted as resulting from a delay of pulses in the counter 
under consideration with respect to the moment of transition of the ionizing particle. 
A pulse delay effect was also observed to occur in similar counters by Crain and 
Griffin (1958). 

In his paper (1958), the author stressed that this delay effect can serve for ex- 
plaining the differences in the results of T}; measurements obtained with the CO, + 
+ CS, counter. 

It is the aim of the present paper to discuss the influence of the delay effect on 
the accuracy of the device used for measuring very weak 14C sources, and to attempt 
an evaluation of Tij from the measurements by Miller and Hawkings et al. taking 
into account the effective time of measurement in their experimental conditions 


as resulting from time distribution data for pulse delay from CO, + CS, counters 


(Mościcki, Renk). 


a) Low level counting 


Let N, and Ly denote the number of ionizing particles per unit time passing 
through the anticoincidence shield and the measuring counter, respectively. Moreover, 
Lo = Bo + 1,9, wherein the index k stands for the penetrating cosmic rays compo- 
nent, and By is the counter background due to radioactive impurities in the counter 
and surroundings. If the a.c. channel gate width is +, then the effective recording 
time 1s; 


x 1 
nope us 


Tea 


and the number of pulses of the a.c. shielding recorded 


No 


N-——ÉÀ—i 
l + Not 


If the usual single channel recording devices are used, the effect of the back- 
ground recorded will consist of the following three components: 


1) Pulses B, not occurring in the dead time of the device 


Bo 


ATTEN z, 


2) Pulses due to particles of the hard cosmic rays component passing in coincidence 


with some pulses No, but falling outside the interval of blocking of the a.c, 
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device as a result of the delay effect. The number of such pulses amounts to 


Log 


ds. 14 Nc 


à (1), (2a) 
wherein ®(r) is the probability of firing the counter after the time r. 

3) Pulses due to particles of the Ly, component passing through the measuring 
counter when the a.c. system is blocked, at a moment 0 < t < + subsequent to 
the front of the a.c. system gate, with delays exceeding t-t, (r denotes the gate 
length). The background due to these pulses will amount to 


Ket No ` Lox è 
0 


The joint background effect arises by summing of the foregoing effects. As proved 
by the results available, and especially by the excellent plateau of this kind of counters, 
the contribution of double and multiple pulses to the total background level is not 
essential. Thus, the total background measured per unit time will be 


B = (1 + N, 1) {By + Loy [Ô (1) + No | D 6— t) dz). (3) 


Assuming a delay distribution law of the form P(t) = ke-^, the following expression 
is obtained for the increase in the number of background pulses (A B) due to the delay 
effect : 


AB (lids La. Roms ead al (4) 
Pune e.g. N, = 10% cpm, r =8.3:x-10% min, 1, = 250 cpm, we have, with 
k = 0.65, 4 = 0.7 x 105 min-!, AB = 2.5 cpm, whereas with A = 0.4 x 105 min-!, 
k= 0.5, we have AB = 7 cpm. 

Of course, an increase in background level occurring simultaneously with 
reducing the time of measurement due to the necessity of applying very long time 
delays of the a.c. shielding channel tends to diminish the accuracy and sensitivity 
of the measurement. A criterion for the accuracy of the measurement may 
be given by C? = S?/B (S denotes the activity of the sample in cpm). In the case of 
a detector for which the blocking time of the device can be neglected, Cj = S2ç/B,. 
For a detector with delayed pulses, in addition to the increase in B according 
to eq. (4), there is a decrease in the number of pulses from the sample investigated, 


go that S = Sv/(1 + No) and * 


À = Fo Lo c DC). (5) 
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The shape of this function for By = 10, Ly = 200, N, = 103, A = 105,k = 1 (graph a) 
and 4 — 5 x 105 min-!, k = 0.5 (graph b) is shown in Fig. 1. 

The possibility of restricting the harmful influence of the pulse delay depends 
upon the use of electronic circuits for blocking the channel of the measuring counter 
for a period v subsequent to the apearance of each pulse of the a. c. shield’). 


2 
(2) 
d a 

D a 

2 a 

15 

1 LA LL AAAA 
D "4 6 8 10 12 14 16 18. DARF 
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Fig. 1 


In this case, the effective time of the measurement is being shorted and is 
satisfactorily approximated by the following expression: 
2 
TA (Bs) 
, € 
Ter = 2 


However, if a circuit of this kind is used, practically all the pulses Z,, appearing with 
a delay lesser than 7 undergo quenching in the a.c. system. 
Here, too, three sources contribute to the background effect. These are: 
1) the background recorded during the effective time of the measurement, 
2) pulses delayed by more than r, 
3) pulses delayed by more than 7 + t (where t < c is the a.c. gate length). 


1) Indeed, such a device was applied with good results by Crane and Griffin. As these authors give 
no details, it may be useful to point out the following possibility. The pulses from the pre-amplifier of 
the a.c. shield trigger a bistable univibrator of appropriately chosen time delay, distributing the pulses 
over two identical channels. The pulses from either channel superpose in the summation circuit, yielding 


a pulse at the output lasting from 7 to 2t according to the degree of superposition of the pulses coming 
from both channels. 
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Neglecting highly improbable cases of more than two pulses occurring within 
the interval 7, the total background is now given by the following expression: 


BY = Bo Tę + Lox (L + No 2)? {8 (9) + No f D(t +0 dą. (6) 
0 


Fig. 1 (graphs c,d) shows the dependence of the ratio C2/C? on 7, with assumptions 
similar to those made previously. 
b) 71 C measurements 

Available 74, 1C measurements include a group of results obtained with 
a CO, + CS, counter working in the G.-M. region. The latter results differ obviously 
from those obtained by other methods?) They are, however, consistent in yielding 
Tij values ranging from 6300 to 6100 years. As the respective authors do not take 
into consideration the pulse delay effect in counters of this type and the related dead 
time, it is of interest to investigate the corrections resulting from this effect relating 
to measurements by the technique considered. T, measurements are characterized 
by a relatively high number of countings, amounting to several thousands cpm. In 
these conditions, the dead time of the counter is quite apparent and affects essen- 
tially the Tj, results. Two of the published available, papers bring material that 
is adapted for checking, namely, those by Miller and Hawkings et al. 

In computing the corrections, the delays distribution was assumed to be of the 


form 
© ()=6 s with a 4s 5x lIütmm.t. (7) 


The value of A=5 x 10% represents a compromise between the value of 
4.2 x 104 obtained by direct measurement and that of 5,1 x 104 resulting from the 
discussion of the data obtained by the author with an M2 counter and taking into 
consideration the dependence of the background on the a.c. channel delay time 
(Mościcki and Renk 1960). 

In Table I, the data published by Miller and the results obtained with the effective 
measurement time correction are assembled. The columns contain: 1— the 
number denoting the counter in Miller's paper, 2 — N — the number of cpm measured 
in that paper, 3 — N,, the number of cpm corrected with respect to the function 
© (t), 4 — the activity of the sample, 5 — f, = No/N, 6 — f — the linearity correction 
factor given by Miller, 7 — f/f. 

Columns 8 and 9 bring the values of the specific activities as obtained by Miller, 
account being taken of the linearity correction factor and the activity in cpm/mol 
with the correction for the dead time of the counter. 

It is to be regretted that Miller does not give the method used in computing 
the linearity correction factor. It seems, however, that his method of determining 
the factor involved systematic error amounting to a coefficient of about 1.037. The 
dispersion of the ratio fo/f in column 7 of Table I is, in fact, small. There is only 


2) For tabellarized results, see e.g. Libby (1952) 
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Table I 
J No à 
Nr N No N,-T | f= W iy, Toll, A Miller | A corr. 
1 2 3 4 5 6 ji 8 9 
18 2090 2151 1861 1.043 1.000 1.043 44.0 46.2 
19 1420 1462 1422 1.030 1.000 1.030 51.0 52.5 
20 2200 2301 2222 1.046 1.005 1.041 47.3 49.1 
23 5230 5841 5682 1117 1.065 1.048 41.3 49.8 
22 7130 8316 8217 1.166 1.138 1.025 46.3 ,47.6 
18 9840 12251 11950 1.245 1.203 1.035 48.8 50.8 
19 9550 11805 11866 1.236 1.295 0,955 (47.4) (45.7) 
47.4 49.3 
(48.8) 


one exception, namely, the last measurement for which the linearity correction factor 


computed by Miller presents an exceedingly high figure. For this reason, this last 


result seems to be less reliable than the remaining ones, notwithstanding the fact 


that the specific activity of the sample as obtained by Miller is in good agreement 


Table II 


A. Set of thick counters of volume difference AV = 15.09 cm? 


V Long counter Short counter Differences 
Volts cps cps corr. cps eps corr. cps cps corr. 
i} 2 4 
2480 37.88 39.68 13.02 13.23 24.86 26.45 
2520 38.48 40.34 13.24 13.45 25.24 26.89 
2560 39.28 41.22 13.54 13.76 25.74 27.46 
2600 39.29 41.23 . 13.84 14.07 25.45 27.16 
2640 39.48 41.44 14.16 14.40 25:32 27.04 
2680 40.11 42.14 14.28 14.53 25.83 27.61 
2720 40.24 42.28 14.84 15.11 25.40 Zi 
25.41 + 0.10 cps 
27:11 = 0.El cns 
B. Set of thin counters AV = 10.79 cm? 
2440 28.00 28.97 9.22 9.32 18.78 19.65 
2480 28.10 29.08 9.08 9.18 18.02 18.90 
2520 28.42 29.42 9.61 9.72 18.81 18.70 
2560 28.61 29.63 9.74 9.85 18.87 19.78 
2600 29.01 30.06 9.84 9.96 19.17 20.10 
2640 29.02 30.07 9.91 10.03 19.11 20.04 


18.76 3:19.53 £0.09 F 


F0.08 
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with the figures for the remaining results. The average of the first six measurements 
yields a figure exceeding the result measured by c. a. 4%, thus leading to a Tj result 
by 4% lesser. By Miller's measurements, Tjj = 6100 years. With the correction, 
Tij = 5864 + 200 years (or 5925 + 250, if seven measurements were taken into 
account). 

High precision measurements of Tip ŚĆ are due to Hawkings et al, who used 
sets of two counters differing in length and working in identical conditions. The device 
eliminated end effects automatically. However, in evaluating the correction for the 
dead time of the quenching system, the foregoing authors didn't take into consideration 
the delay between the moment of transition of the ionizing particle and the *discharg- 
ing” ofthe counter. The results obtained by them, and the results corrected, are 
given in Table II (detailed results of one measurement) and Table III (ultimate 
results of three measurements). 


Table III 
Measure- Activity exp. Activity corr. 
eR Set of counters (4) (4,) As A, 
I Thick 1.684 1.798 1.068 
Thin 1.738 1.842 1.060 
II Thick 0.868 0.890 1.025 

Thin 0.860 0.889 1.034 

II Thick 1.298 1.366 1.052 
Thin 1-254 1.300 1.034 

1.046 


It is seen from Table II that, when neglecting the effect of delayed pulses, the 
authors obtained figures for the activity of the samples that were approximately 
5% too low. This leads to the corrected value of T,4C as given by the foregoing 
authors of 6320 + 200 to about 6040 + — 300 years. The preliminary character 
of the ®(z)-function can produce a difference in correction-numbers amounting to 
10—20%. Nevertheless it seems that the pulse-delay effect in CO, + CS, counters 
can quantitatively explain the discrepances in T}; measurements. 
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ÉQUATIONS DU POTENTIEL ET DE LA TEMPÉRATURE POUR 
CONDUCTEURS DANS UN CHAMP MAGNÉTIQUE 


Pan H. SrACHOWIAK 
Institut de Physique de l'Académie Polonaise des Sciences, Wroclaw 
(Recu le 19 Juillet, 1960) 


En imposant les conditions stationnaires aux équations de la thermodynamique 
des phénoménes irréversibles concernant les effets thermo-galvanomagnétiques on obtient 
les équations du potentiel et de la température. Celles-ci donnent la solution compléte 
du probléme pour toutes les conditions aux limites, contrairement au formalisme jusqu'ici 
en usage. On a prouvé qu'il suffit de donner la valeur de la température à la surface pour 
obtenir une solution unique dans le cas où le conducteur est en isolation électrique. On 
a donné une méthode de trouver cette solution. On a montré que dans un tel conducteur 


peut apparaitre un courant électrique stationnaire. 


1. Introduction 


L'article présent est le développement de deux articles publiés dernierement 
(Stachowiak 1960). 

Avant l'origine de la thermodynamique des phénoménes irréversibles, la théorie 
des effets thermo- et galvanomagnétiques était basée sur une méthode inexacte dite 
quasithermodynamique (voir de Groot 1952, Callen 1949). 

. Callen (1948) a élaboré la théorie de ces effets fondée sur la thermodynamique 
des phénoménes irréversibles. Les équations phénoménologiques obtenues expriment 
le courant électrique et le courant d'entropie comme formes linéaires des composantes 
du gradient du potentiel électrochimique et de la température. 

Mazur et Prigogine (1951) ont obtenu des équations un peu différentes, ou le 
gradient du potentiel électrique et le courant de chaleur sont exprimés par le courant 
électrique et le gradient de température. 

Ce probléme a été aussi considéré par Fieschi et autres (1954) et de facon rela- 
tiviste par Kluitenberg et de Groot (1953, 1954, 1955). 

Dans l'article présent nous emploierons le formalisme de Callen généralisé sur 
trois dimensions (voir Landau et Lifshits 1957) ; celui-ci offre le passage le plus commo- 
de aux équations du potentiel et de la température. 

(67) 
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Les équations phénoménologiques, dans notre cas au nombre de 6, relient entre 
elles 12 grandeurs (6 composantes des courants et 6 composantes des gradients) et 
par conséquent elles ne peuvent donner la résolution compléte du probléme. Nous 
devons connaitre 6 grandeurs parmi 12 pour fixer les 6 autres. Le but de cet article 
est d'obtenir les équations du potentiel et de la température qui donnent la résolu- 
tion complete du probléme pour toutes les conditions aux limites. Le formalisme des 
équations phénoménologiques employé jusqu'ici par les auteurs ne le permettait 
pas. Nous considérons uniquement les effets linéaires par rapport au champ magnétique. 
Le conducteur est supposé non ferromagnétique et les coefficients d'Onsager dépen- 
dants de la température seule. Le cas d'un conducteur en isolation électrique va nous 
intéresser particulierement. 

En $2 nous obtenons les équations du potentiel et de la température et nous 
examinons les conditions aux limites. En $ 3 une méthode de résoudre ces équations 
est élaborée et la preuve de l'unicité des solutions est donnée. Le $ 4 est consacré 
au probléme de l'existence d'un courant électrique stationnaire. 


2. Les équations du potentiel et de la température 


Nous allons nous servir des courants et des affinités employés dans la mono- 
graphie de de Groot (1952). Alors les équations phénoménologiques à trois dimensions 
pour un conducteur placé dans un champ magnétique ont l'aspect (comp. Landau 


et Lifshits 1957) 
-j= LVu +h VT + hs [Fa B] + [V T, B]. (la) 
— Š= ha Vu + LPT + ha [Vm B] + la [V T, B], (1b) 


où J est la densité du courant électrique, 5 — la densité du courant d'entropie, B — le 
vecteur de l'induction magnétique; 4, ły, Ly, lig, La et lag sont les coefficients d’Onsa- 
ger: particulièrement ły est le coefficient de conductibilité électrique, LT celui de 
conductibilitć thermique; m est le potentiel électrochimique. 


B = B, t M; (2) 


où u, est le potentiel électrique, u, — le potentiel chimique. 
Les équations (1) négligent les effets dépendants des puissances supérieures du 


champ magnétique. C'est une approximation juste pour les champs magnétiques 
faibles. 


L'application des conditions stationnaires aux équations (1) donne les équations 
du potentiel et de la température. Ces conditions sont les suivantes (voir Callen 1948): 


* 


p-J=0, (3a) 
p 


© 
I 


T | (3b) 


69 


= 
ou Q est la densité du courant de chaleur 


OE TS. (4) 


Les équations (1) donnent 
—p-F=h du + AT—lspulp. B— lapTip. B] HL pT pu + 
+ p (p T? + hs B[p T, pal, 
—r-0—ül TAu ETAT IL Tpulp, B] — 4 TpT[p, B] + 
+ (al) y Typ + (LTY (p TP + Gal) b [y T. pu], 


oü l'apostrophe indique la dérivée par rapport à la température. 


> > 
Nous pouvons négliger les thermes contenant rot B, car B est de principe un 
champ extérieur, done son rotationnel doit étre zéro à l'intérieur du conducteur 
(nous ne nous occupons pas de l'influence que peut avoir le courant électrique présent 


dans le conducteur sur B). 
Si le courant électrique est proportionnel à B, ce qui a lieu en effet, comme 
nous allons voir dans le paragraphe suivant, lorsque le conducteur est isolé au point 


de vue électrique, on obtient de l'équation (la) 


hyu + hapT-B. 
Donc 


We y T] — B. 
Négligeant les thermes bilinćaires en B nous obtenons 
ax Mi l Au + la AT + hp Tp ui + ha (p TY, (5) 
—p:Q- d TAa STAT + (là TY pT yu + 0 T (p T). 
Les conditions (3) nous donnent donc les équations suivantes: 
h Au + he AT + Lp Tu + ls (p TP = 0, (6a) 
la TAu + STAT +h, (Pro? + tha + Ga T) Tye + (TY (p T? = 0. (6b) 


Les équations (6) ne contiennent pas B. Donc dans l'approximation que nous sulvons 
les équations du potentiel et de la température ne dépendent pas du champ magnéti- 
que, si le conducteur est en isolation électrique. Les conditions aux limites pour un 
tel conducteur ont la forme 


T (c) = f (o), (7a) 
EAM. (7b) 
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où o indique un point de la surface du conducteur et J, la composante normale du 


= PY. — b. . 
courant J à la surface. Évidemment les conditions (7) ne sont pas les seules possibles. 


Nous les choisissons parce qu'elles ont une interprétation physique simple. 

Occupons nous de la condition (7b). Choisissons dans un point donné de la surface 
un systéme de coordonnées cartésiennes 7, t, v, tel que l'axe n soit normal à la surface 
et les axes t et v tangents. Alors la condition (7b) prendra la forme 


0 Ou Ou 
I a + hs By D —hB 2, 


T oT 
=the ot lu Be zę = ee) (8) 


ou autrement 


y l Ou 
2 2 2 | 2 y, A 1 | 
lis Be ou L = hs B, ar) 
TR BB BS ot | [2+ B,(B? + BS[ 8v 


l aT 
mc n 2 2 21114 12. 
= — [If + lia (Be + Bo) m TR (B? IBS == Jie 


lig Bo zur! — l B: A (9) 
[ie + It (B; + Be ))^ ot T IRSE DEN (B? + Bi] dv 


Nous avons donc une égalité entre les dérivées suivant diverses directions. Négligeant 


les thermes carrés en B nous pouvons écrire 


PR a di (10) 
os] L Is,” 
qe 
où — et — indiquent les dérivées suivant les directions définies par 


95, OSs 


> l m > 
sa mg Bb uno, (11) 
12 12 


> 
n, t, v sont les vecteurs unitaires répondant aux axes m, t, v. 

Nous voyons que le champ magnétique a une influence sur la forme des condi- 
tions aux limites, Si B disparaît nous avons un probléme de Neumann, en cas contraire 
nous avons une condition à dérivée oblique, telle que l'angle entre la direction oblique 
et la direction normale dépend du champ magnétique, 
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Nous aurions pu imposer des conditions aux limites indépendantes du champ 
magnétique. Mais dans ce cas le changement du champ magnétique changerait leur 


sens physique, car B entre dans les équations (1). 


3. Résolution des équations du potentiel et de la température et problème d'unicité 
A) en absence du champ magnétique 


Avant de nous occuper des équations (6) en présence du champ magnétique, 
nous devons prouver que la condition sous laquelle elles ont été obtenues, c'est-à-dire 


que J est proportionnel à B, est satisfaite. J(x, y, z) est une fonction du champ magnéti- 


que. Introduisons une unité d'induction By(x, y, z) telle que 


B(x, y, z) = B Bo, y, 2), (12) 
ou ß est la mesure du champ magnétique. 


Si le champ B est faible nous pouvons développer J en série de Taylor: 


= © + à 6 B) B + 36. Bo) P + ... (13) 


> 


ją et ją sont évidemment des fonctions du lieu et des fonctionnelles de By. Le premier 
therme à droite est indépendant du champ magnétique, le second est linéaire en B, 
le troisiéme bilinéaire. Si nous prouvons que J, = 0, J peut étre en effet considéré 


comme proportionnel à B dans notre approximation. 
S'il ny a pas de champ magnétique, les relations phénoménologiques ont l'aspect 


+ 


—J = LV wu + ls V T, (14) 
OL. True IS TV T. (15) 


En inroduisant (14) et (15) dans (3) nous obtenons les équations (6). Les conditions 


aux limites auront la forme 


+ 


T(o) = f(o), (16) 


Les solutions de ce probleme sont des systómes de fonctions u(t), T(r). Nous allons 
démontrer qu'une fonction T(r) donnée ne peut entrer que dans un seul systéme 
u(r), T(r) constituant une solution, c'est-à-dire que si nous connaissons T(r) qui 
est l'une des fonctions dont se compose la solution, la seconde fonction p (r) est 
unique, à part une constante additive. En effet, si T(r) est une fonction donnée, 


l'équation (6a) est du type 


Au Y fo o0 (18) 
j=1 
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avec la condition aux limites 


25 h (o) (19) 


Ce probléme n'a qu'une seule solution, mettant à part une constante additive (voir 


Miranda 1955). 


Il est facile de voir que w remplissant la relation 
l 
yu—-—-pT (20) 


est une solution particulière de (6a). Puisque ly, 1, dépendent uniquement de la 
température, la part droite de (20) peut étre écrite sous la forme d'un gradient. La 
relation (20) définit w univoquement, à part une constante additive. y satisfaisant 
à (20) satisfait aussi la condition (17) et constitue par conséquent la seule solution 
du probléme (la constante additive n'a pas d'influence sur la relation (20)). En introdui- 


sant (20) dans (14) et (15) nous constatons que JE 0, et 


SG UI RT (21) 


(L l — R) T 
lą 


C(T)= (22) 


Introduisant (21) dans (3b) nous obtenons 


p : (C(T) VT) = 0 (23) 


avec la condition aux limites (16). Ce probléme n'a qu'une seule solution, car nous 
pouvons passer de la variable T' à la variable O définie par l'intégrale (voir Carslaw, 


Jaeger 1959) 


TQ) 
l 
© (r) = —~ C (T) aT. 
T 
ou 7 a une valeur queleonque mais fixée. Introduisant (24) dans (23) nous obtenons 
l'équation de Laplace 


AO — 0 (25) 
avec la eondition aux limites 
f(o) 
il 
© (o) = C) fem als (26) 


Ce probléme bien sür a une solution unique. Le passage de © à T est univoque, si 
Ja Y . A . . . Z 
le coefficient C(T) a toujours le méme signe. Le deuxième principe de la thermodyna- 
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mique appliqué à (21) fait que cette condition est remplie (la direction du courant 
de chaleur est inverse à celle du gradient de température). 

Nous voyons done 

1) que le probléme (6), (7) en absence du champ magnétique n'a qu'une seule 
solution. 

2) qu'il n’y a pas de courant électrique. 


B) en présence du champ magnétique. 


Nous avons le systéme d'équations (6) et les conditions aux limites (7a) et (10). 
La relation (20) satisfait toujours à l'équation (6a), mais elle ne remplit plus la condi- 
tion aux limites, car s, et s, dans (10) n'indiquent pas la méme direction. Nous divi- 
sons # en deux parties selon la formule 


yu=— pT+ py. (27) 

Introduisant (27) dans les ćquations (6) nous obtenons 
L 4o + LV TV 9 — 0, (28) 
V (CT) VT) + V: (ls Tg) — hayTyo + Lmg) —0 (29) 


avec les conditions aux limites 


T(o) = fto). (7a) 


9p  L 9f 9f 30 
Da h (2.2 Bia. i EU, 


où 
ré ls hu (31) 
h 
Apercevons que la part droite des conditions aux limites ne contient maintenant que 


des fonctions connues. 
Nous passons à la variable © à l’aide de la transformation (24), ce qui permet 


d'écrire 


dasha 
C(z)4 8 = —4 35 A T(0())yO yo (33) 


avec les conditions aux limites (26) et (30). Le point indique la dérivée par 


rapport à O. 
La relation (27) prend la forme 


Le CG) 
Fu =— T co) V9 * vv. (34) 
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O (r) et p (r) peuvent être développés en série de Taylor selon B, ainsi que J 
l 


0 (7) = Oy () + BY (1). (35) 
p) = B qi (1. (36) 


Les thermes non linéaires en B ont été négligés. Introduisant (35) et (36) dans (32) 
et (33) nous obtenons dans l'approximation linéaire en B 


AGO, —U, (37) 
Ag + lp Or 9 — 0, (38) 
d l 

CFO =—h = T (8%) 7 Oo yg. (39) 

dO, V d, 

ou ©, est défini selon la formule 

0 = 6, + 6,. (40) 

Les conditions aux limites sont les suivantes: 

fio) 


O, (c) = | V JE (41) 


C Le 
dp OL of of 
9n |, (s. at e dv}? ies 
©, (c) = 0. (43) 


Nous résolvons d'abord l'équation (37) sous la condition (41). Ce probléme 
n'a qu'une seule solution. Puis nous introduisons 6% ainsi obtenu dans l'équation 
(38) et nous cherchons la solution de cette dernière sous la condition (42). Diverses 
solutions d'un tel probléme ne peuvent différer que d'une constante additive (voir 
Miranda 1955). Ensuite nous introduisons @ et 6% dans l'équation (39). (39) et (43) 
donnent une solution unique, la constante additive dans g ne joue évidemment aucun 


róle. 


4. Le rotationnel du courant électrique 


+ 


Puisque la divergence de J disparaît et qu'il n'y a pas d'afflux de l'extérieur, 
sil y a un courant électrique, il doit y avoir des lieux où rot JÆ 0 (à moins que le 
conducteur ne soit pas simplement connexe) et réciproquement, si rot J n'est pas 
zéro dans tout le conducteur, il doit y avoir un courant électrique. Calculons donc 
rot J à partir de l'équation (la), en nous servant de (27): 

> " je " -> 
- Ir A=4 lT pol — L [yT,[yT,B]] + 
HEBAT- L(Bp) pT+L(pTy)B. (44) 


Si nous négligeons la dépendance thermique des coefficients et appliquons un champ 


magnétique homogéne 


Br MT 
(45) 


= 

eee ee T 

JE on 

A £ C E 1 | WA pres ? A * D 

où 6,7, 6 sont des coordonnées cartésiennes telles que l'axe ¢ est parallèle à B. 

Nous voyons que, si les secondes dérivées de la température ne disparaissent pas; 

méme dans un conducteur isolé au point de vue électrique peu apparaître un courant 
électrique stationnaire. 


Conclusions 


L'application. des conditions stationnaires aux relations phénoménologiques 
concernant les effets thermo-galvanomagnétiques méne aux équations du potentiel 
et de la température. Ces équations, dans l'approximation linéaire ne contiennent pas 
le vecteur de l'induction magnétique. Le champ magnétique a par contre une influence 
sur la forme des conditions aux limites. Parmi les conditions possibles nous choisissons 
celles qui correspondent à l'isolement électrique du conducteur, avec une répartition 
fixée de la température à la surface. 

Nous avons obtenu une méthode de résolution de ce probléme dans l'approxima- 
tion linéaire en B. Nous avons démontré que dans cette approximation le probléme 
a une solution unique, mettant à part une constante additive dans le potentiel. Les 
approximations supérieures de la solution seraient dépourvues de sens, car les rela- 
tions originelles ne sont exactes que dans l'approximation linéaire. 

Le calcul du rotationnel du courant électrique montre que méme dans notre 
cas d'isolement électrique, un courant électrique peut exister, si les secondes dérivées 
de la température ne sont pas zéro. 

La résolution des équations du potentiel et de la température et un examen plus 
minutieux des conditions d'existence d'un courant électrique stationnaire dans les 
cas les plus simples auront lieu dans un article prochain. 

Je voudrais remercier M. le Professeur R. S. Ingarden pour lintéressement 
bienveillant et continu qu'il a manifesté à ce travail, M. le Professeur W. Rubinowicz 
qui a attiré mon attention sur la question d'unicité et M. le Professeur J. Łopuszański 
pour nos discussions nombreuses. Je remercie aussi M. le Professeur B. Makiej et 
M. R. Freud qui m'ont tenu au courant de leurs expériences et grâce auxquels je me 


suis intéressé à ce probléme. 
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PROTON SPIN-LATTICE RELAXATION TIMES IN VERY DI- 
LUTE AQUEOUS SOLUTIONS OF Ni** AND Mn** IONS 


By J. W. HENNEL, K. KRYNICKI, T. WALUGA, AND G. ZAPALSKI 
Institute of Nuclear Physics, Cracow 
(Received October 2, 1959) 


Proton magnetic spin-lattice relaxation time T, in air-free aqueous solutions in the 
Ni** concentration range between 2x 1015 and 2X 1018 ions/cm? was measured. For Ni** 
and Mn** ions the dependence of l/T, on concentration troughout the temperature 
range 10 — 90°C was shown to be linear. A simple method was developed for 
obtaining air-free samples of aqueous solutions. 


1. Introduction 


The addition of paramagnetic ions to water decreases the proton spin-lattice and 
spin-spin relaxation times, 7, and 74. This effect was first interpreted by Bloembergen 
et al. (1948) in terms of an ion-proton magnetic dipol-dipol interaction, which predicted 
a proportional dependence of 1/7, on the concentration: 


1/T, = ke, (1) 


where c is the concentration and k is a function of the temperature only. This depend- 
ence is found to be correct for not too weak solutions. 

In the case of very dilute solutions the interaction between protons must also 
be taken into account. As 1/7, is proportional to the transition probability and because 
of the additivity of transition probabilities resulting from the two interactions, the 


expression for 1/7, is of the following form: 
Ulli 2 he, (2) 


T,,, represents the proton spin-lattice relaxation time in pure water. 

Recently new theories of nuclear relaxation processes in paramagnetic ions 
solutions were developed (Bloembergen 1957, Hausser and Laukien 1959, Bernheim 
et al. 1959). These theories, successful in explaining the inequality of T, and T, in 


(72) 
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some ion solutions, also take into account scalar coupling between the nuclear spin 
of the proton and the electron spin of the paramagnetic ion. With respect to these 
theories, eqs. (1) and (2) still remain valid. 

In the range of medium and strong concentrations of paramagnetic ions, Tj 
measurements were made by a number of authors. The most extensive investigations 
of the concentration dependence of T, (and T3) were carried out by Laukien and Schlüter 
(1956), Nolle and Morgan (1957) and Laukien and Noack (1958). The results obtained 
prove the validity of eq. (1). 

In the low concentration range, where T, > 0.1 7,,,, only few results for air-free 
solutions are known. Giulotto (1956) carried out measurements of 7, in water and 
some organic liquids containing the paramagnetic gases O, and NO in solution. 
His results corroborate the validity of eq. (2) in this case. 

In the present paper measurements of T} in the above mentioned low concentra- 
tion range of paramagnetic ions Nitt and Mn** are presented. These ions were 
chosen because of their different T,/T, ratio behaviour and absence of hydrolysis 

effects. 


2. Preparation of Samples 


As the atmospheric oxygen dissolved in water affects 7, considerably, it was 
necessary to develop a method for preparing air-free solutions. This was achieved 
by boiling at length in the simple glass apparatus shown in Fig. 1. 

A known amount of paramagnetic salt in aqueous solution was introduced into 
the glass bulb (a), which was then filled up with carefully purified water. The inlet 
tube was sealed off at point (b). By heating the bulb, the solution was evaporated to 


Fig. 1. Glass apparatus for preparing air-free solutions. 
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approximately 1/4 of its initial volume. The steam carrying away the air in solution 
escaped through tube (c), the end of which was immersed in water. On switching off 
the heater, the water level in tube (c) began to rise and at this moment the apparatus 
was sealed off at (d) yielding a closed system containing the air-free solution. 

The arrangement together with the sealed off tubings was weighed to determine 
the concentration of the solution. Subsequently the solution was poured into the small 
bulb (e) constituting the future sample, which was then sealed off at (f). 

For preparing the solutions, A.R. grade NiSO, : 7H,0 and MnCl, : 4H,O were 
used. Water was purified by the following method: ordinary laboratory distilled water 
was treated with 25 x 10-? g KMnO,, 1 g AIK(SO,), - H,O and 0.7 g NaHPO, + 12H40 
per litre and then distilled twice: 

High purity of the glass apparatus used was most essential. 


3. Experimental Procedure for T, Measurements 


A method described elsewhere (Hennel and Hrynkiewicz, 1958) was used for T, 
measurements at 28 Mc/sec. 

For detecting the nuclear magnetic resonance, the high frequency bridge described 
by Grivet et al. (1951) was used. The bridge was mounted in a box made of a brass 
block. The entire bridge device containing the sample was kept at the constant temper- 
ature required. This was carried out by letting water flow from a thermostat through 
copper tubings soldered within the block. 

The temperature was measured with a glass-mercury thermometer mounted in 
the outlet tube of the block. The readings of this thermometer were calibrated by 
comparison with an exact resistance thermometer situated in place of the sample. 
The corrections obtained amounted to (0 + 0.1)?C for 20°C and (1.5 + 0.3)°C for 
90°C, and were reproducible. 


4. Results and Discussion 


The resulting values of T, for 20?C are shown in Fig. 2 versus the Ni** ions con- 
centration. The theoretical curve drawn according to eq. (2) assuming! Tie = 3.1 sec 
and k = 1.02 x 10-18 cm?/ion. sec is also shown. This value of k at 20°C was calculated 
from the data obtained for the much greater Ni** concentration of 5.5 x 10% ions/cm? 
by Hausser and Laukien (1958) at 26.5 Mc/sec. 

All T, measurements in the present investigation were accurate to 3%. 

The experimental results in Fig. 2 support eq. (2) and show that k remains con- 
stant throughout the concentration range of 2.38 x 1016 to 5.5 x 10% ions/em?. 
This can be stated to an accuracy of about 10%, except for the lowest concentration 
range of about 3 x 1016 ions/cm?, where the accuracy is 30%. Taking into account 
the measurements of Laukien and Schlüter (1956), the constancy of k at room temper- 


1) The interpolated value of 74 for water at 20°C is from Hennel et al. 1958. Other values are 3.1 sec 
(Giulotto et al. 1957), and 2.92 sec (Simpson and Carr 1958). 
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ature throughout the Ni^* concentration range of 2.38 x 1015 to 2.3 X 10% ions/cm? 
has been proved. 

To investigate the validity of eq. (2) in the case of Mn** ions, two samples contain- 
ing 1.19 x 1016 and 1.33 x 1016 Mn** ions/cm? were prepared. The measured 77 
values at 20°C were 1.97 sec and 1.92 sec respectively, and the calculated k values — 


10 10 10° ions/cm? 10° 


Fig. 2. Proton relaxation time T, versus the concentration of Ni** ions in oxygen-free water at 20°C 
o experimental, — theoretical. 


1.55 x 10-17 and 1.49 x 10-17 em? per ion and sec. This is in good agreement with 
the value k = 1.46 x 10-17 obtained at the much higher concentration of 3 x 1018 
ions/cm? by Bernheim et al. (1959, Fig. 1) at the same frequency and temperature. 

To investigate the validity of eq. (2) in the temperature range between 10°C 
and 90°C, three of the afore mentioned samples (two containing Nitt ions and one 
containing Mn** ions) were used for 7, measurements at different temperatures. The 


results are compared with 7, values found at much higher concentrations by calculating 
the expression: 


Cas P a el 35) lin" fp 


where 7°, is the spin-lattice relaxation time at high concentration c. This expression, 
according to eq. (2), should be independent of the temperature and c4, equal to 
the concentration of the sample. This is shown in Table I and proves the validity of 


eq. (2) throughout the temperature range 10—90°C. To obtain c, we used T, 
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values interpolated from our experimental results shown in Fig. 3. The Ty, values 
are taken from Hausser and Laukien (1959) at c, = 5.5 X 10% Nitt ions/em? and 
from Bernheim et al. (1959) at c, = 3 x 1018 Mn++ ions/cm. 

Fig. 3 shows the temperature dependence of the function A — T, n/T (wherein 7 
is the viscosity of water, and T — the absolute temperature) as found experimentally 


Tm 
T 
sec. porse 


10 20 30 40 50 60 70 80 90 96 


Fig. 3. Temperature dependence of the function T,5/T for very dilute aqueous solutions of para- 
magnetic ions: 


(a) 1.33 x 1016 Mn** ions/cm, (b) 1.53 x 1017 Ni ions/cm?, (c) 6.89 x 1017 Ni** ions/cm?, and 
(d) pure water (Hennel et al. 1958), for comparison. 


for the three solutions of paramagnetic ions. This function should be constant according 
to Bloembergen et al. (1948), but, as a matter of fact, it decreases towards higher tem- 
peratures even in pure water (Giulotto 1957, Hennel et al. > Simpson and Carr 
1958). It is seen from Fig. 3 that a small amount of paramagnetic ions causes a con- 
siderable rise in the slope of the 4 curve. This can be explained by the temperature 
dependence of k for the ion under consideration. | | 

The foregoing fact can have some influence on the unexplained behaviour of A 
in pure water, as even in very carefully purified water the presence of some slight 


paramagnetic contamination is possible. 
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Table I. 


Concentration coace = (1/T1—1/Tiw) Tin’ cj as calculated from Tj at various temperatures (°C). Inter- 


posten values of Tą are given in Table II. The concentration is in 10" ions/cm? units. 


|ESES="IN “Calculated concentration MEA t 
Ion" concentration 10° | 15° | 20° | 40° | 50° | 60° | 80° | 90° h 
| | 
Nit+ 153 17.1. | 17.12.17.0-232.1.16.97016.97171622116.2221.521816:6 55000 
Nit+ 68.9 62.0 | 62.7 |63.5+6 64.6 .| 64.6 | 64.7 |65.0+4 63.8 55000 
Mn++ eae 1.36 | 1.38 | 1.404-0.3| 1.46 | 1.42 | 1.43 | 1.444-02 1.48 | 300 


Table II. 


Interpolated values of the spin-lattice relaxation time T}, for oxygen-free water at various temperatures 


(Hennel et al. 1958). 


ssij o] 10 | Bo] 0 | 90] 40 | 0 | 0 M TEL 
pe | um 235 | 27 | 81 Mad | 50 | 66 | 2 4 Sun 
(sec) | 
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SUPPLEMENTARY NOTE TO THE PAPER: 
MOLECULAR THEORY OF LIGHT SCATTERING BY MULTI. 
-COMPONENT SYSTEMS 


By 5. KiELICH 
Institute of Physies, Polish Academy of Sciences, Poznań 


(Received October 17, 1960) 


The present note is to supplement the paper on a *Molecular Theory of Light 
Scattering by Multi-Component Systems" (Kielich, 1960), and brings a discussion 
of the scattered light intensity for the case of a system whose various components 
consist of quadrupolar, anisotropically polarizable and hyperpolarizable molecules. 

The light intensity scattered at the angle Ÿ is given by the general formula (see, 
Kielich, 1960) 


824 I 


I) = 45 45 A^ R2 Ri 


(5( (1 + cos? 9) F(s) + (13 + cos? 0) Fanis (s)} (1) 
where I, and A denote the intensity and wave-length of the incident light, respectively 
and R, is the distance of the point of observation from the centre of the scaiter- 
ing volume V. The factors F; (s) and F,,,, (s) characterize the molecular mechanism 
of isotropic and anisotropie light scattering, respectively. For a multi-component 
system these are of the form 


xiN xjN (Psi) (4,3) | * t 
Omg Im, —i8* rbd 
ee > (n tee > 2 IEP” E: ds. ) » 
ij p—l g=1 

“N xN (pst) (43) V * \ 

On Mg Om; sn) 
A (GR 5 D ( (30, Óps = ap Óy6) > = (oa £ 5 

Ly 


(3) 
where x; is the mole fraction of the i-th component, m? — the &-component 


of the Em dipole moment induced by the electric field E, of the incident light 
wave in the p-th molecule of species i, r(?9 —thes vector connecting the centres of the 
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p-th and q-th molecules of the i-th and j-th species, and, for Rayleigh scattering, 
|s| = (4x/2) sin (8/2). The brackets Jin eqs. (2) and (3) denote the statistical 
average in the absence of the external electric field (E = 0) and 9,5 is the substitution 
tensor (unity if a = f and zero when «Æ f). 

The total differential polarizability tensor of the molecule immersed in the medium 
is given by the following expansion (for E — 0): 


PC ; pod lOD pod) pi) 
(one) [ago + aa > à Lg ceo à 
JE?" /E=0 


1 (Psi) | 
+ = Be Fa +. ‘| DES: Te 


4 
3 DEV fi =, ? ( ) 


with qus denoting the polarizability tensor of the p-th molecule of species i, and 
Be vB, —its hyperpolarizability tensors. The tensor BY, accounts for the 
additional polarization of the p-th molecule of species i as idee by the field grad- 
ient F@®. The tensors &,5; Boys Vagyo and B,5.,, have been discussed by Buckingh- 
am (1959). FO) and F& are the a-component of the molecular field and 
«p-component b the Gad, eens at the centre of molecule p of species i due to the 
electric charge distributions of all the others, in the presence of the electric field E, . 

Substituting the expansion of eq. (4) in eqs. (2) and (3), the following expres- 
sions are obtained: 


x,N x; N 
5) = XO Ż («9 al AD „+ HD, cos s + 995 (5) 
iJ xu i 
xN X; N | 
Pas © = 4363) 2 (0 a+ AB 4 
= i qm 
+ He Base ape io —0 5 0,5) cos S * re > (6) 


Here, X225? denotes the cosine of the angle subtended by the œ and y axes of the 


molecular systems of reference rigidly attached to the p-th and q-th molecules of 
species 7 and j respectively. The tensors 


(i) (i) G) 
Aa G) 0) IF G) _G) OF, (i) (D OF, OF, 7 
aß:yd = Cap Kye G) + ys Sas —— (i) F ge Km = @ Gg»? ( ) 
DE, DES SEPOEŃ 
Ha RARE? GRAB. aS O ald + 
) 
3 Tio Bow, ru +1 e Be. :en Ks vs (8) 


determine the effect of the anisotropy of the molecular field and that of the hyper- 
polarizability of the molecules on light scattering in a dense medium. 
With respect to Onsager's (1936) model, the anisotropy of the molecular field 
is given by 
OF, 


0) 
apo. Ye Vas (9) 
"è i 


and the molecular field F;? and field gradient F! 2 at E = 0 are 


pm yË We p go 9% TB 


Qt) edis (10) 
; AW 1 ON) 
FY ciuem ACTE 
aß Org fa m ? (11) 


wherein u® is the &-component of the permanent electric dipole moment of an 

isolated polar molecule of species i and GC? —its electric quadrupole moment. The 

parameters y®, zę and £O are given by: 

2(n? — 1) {(n® + 2) (n? — 1) 4? — (n? + 2) (n? — 1} 

(n? + 2) (2n? + 1) (n? + 2) — 2 (n? — 1) (à — 1) 44)” 
2 (e — 1) (n? — 1) 

(e + 1) (n? + 2) — 2 (e — 1) (a? — 1) 4e” 


i 2e +1 mois 
(1) i a) 
zn 1 
> e eee Ka i) 


where e and n are the electric permittivity and refractive index of the medium, respec- 
tively, and 4? = ala), with a(? denoting the polarizability component in the 
direction of an a-principal axis of an isolated molecule of species i, and 4 — its 
mean. polarizability. 


prom (12) 


xe- (13) 


For non-dipolar molecules (u, = Pepy = 0) possessing the axial symmetry, 


apy 
the following tensor elements only are non-vanishing: 4,5, Og Yapy and Bygiys 


(the axis of symmetry being parallel to 03, see Buckingham 1959): 


Mi emp EA Tue Ya Nr) 5 3Y 1122 = "fo. (ABE == ERES. Alle 


&33 = A), Y3333 — MP 05, = — 20, = —20» = ©, (15) 
Bąą:33 = — Bag = — 2B53:22 = 2 (Biu + Pis) = B 


where © is the quadrupole moment of the axially symmetric molecule and B—its 


quadrupole polarizability. 
With respect to (15) and for Az» ri. 


anisotropic light scattering assume the form: 


the molecular factors of isotropic and 
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x © Ó,, + T y | Geos Ol) E (t) — 7 a, (17) 


86 


which holds for quadrupole molecules having the axial symmetry. Here, ©,, is the 
angle between the axes of symmetry of the molecules of species i and j, r respecte 
g;;(r) — the distribution function and Q*dr = dr,, do; do,, where Q = : (doy. The 
quantities 4@ and Agi), in the foregoing expressions account for the effect of the 


anisotropy of the molecular field on isotropic and anisotropic light scattering, respec- 
tively, and are given by 


i) 1 (9G) gi) L943)040Y 1 1190) of) 1 930 40 
AD =} (AP yO + 240 (9) + à (AT vi? + 24? vf) + 
O0 1. 930 9D) (9G) KG. 4. 930 GO TS 
+3 (AP vip + 22? v3) (AP vi + 22 vi). (18) 
HX B) any eto 16) „A _ 26) „(9 
67 2 A V] zo Wl ER Ap WM NAA 
GA ct (1) 40) (2) 4) 
Aji s AY À AS 
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ME E Vi AL Vl | (A vil TET E (19) 
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HG? and HÈ? account for the additional contribution to isotropic and anisotropie 
light scattering due to higher order effects such as hyperpolarizability of the molecules 
and polarizability of molecular quadrupoles. For these quantities we have: 


n ES) a) VO 10 86 : o 70 FO Ot N 
Hy = 4a is ab wro UD Mis Tape rs (20) 
2 A) A) 
7 ee 4 0. 3 
Ho? = n fan MN "T yt, pia Ya 9 i) 605 ds 
anis anis ns Te n 
Lead qi) Da(i) Am „ad E 
1, By6% ] By, O; 
ai} ai) ad baci) 
Here 
(i) (i) 
: ai —® 
dzy LE and a= 3 (af + 2a) (22) 
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describe the anisotropy of polarizability of the isolated molecule of species “i” and 
its mean polarizability, respectively, whereas 


do- it By = 4y D) 
2174) 


| ra= Omas (3 


are the anisotropy of hyperpolarizability of the molecule of species i and its mean 
hyperpolarizability. 


The parameters 745, Egy xÈ and x, introduced in eqs. (20) and (21) are 
defined as follows 


nen wz d [26x 1l. [nt — IN" gr ro 
n= x (Fe) (3) (xt +2211)» (24) 


2e - 1| [n -1 O oli) | 
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A = (15yg “5 2 GP + 40) Fa — FD} (y FP. (26) 


(i) 1 1 7 m2 = 
anis — {2176 ô y(t) Fis A (350 — 2) (FA — Fb Uva Ów) DE (27) 


T | M ns = ; a 
wherein p mean square molecular fields F4, F; and F“ are given, for axially 
symmetric quadrupolar molecules, by: 


T LC 
FA = FR [x 20 pu n?--2 AE ads (28) 
=) 9 Jet - RET. *[s A O?. 
2, — 1 (i) (i) ; 
FR 5 E = EE Al als , (29) 
(i) 
Sá set, = 0? 
F} = 2F + Fis = 93a? - (30) 
(i) 


In a former paper (Kielich 1960), when passing from the formulas (5.2) and (5.3) 
to eqs. (5.13) and (5.14), the simplifying assumption of xf = x, = 1 had been 
made. From the foregoing eqs. (26) and (27), this condition is seen to be satisfied 


strictly only if the molecular field possesses the spherical symmetry, i. e. if p = 


= F2, = p z prt Assuming also yum ==} à (yf? + yf) (see, Buckingham and 
Stephen 1957), the quantities 0,4 and y defined above by eqs. (23) reduce to 
o b) 
= 767 sO 2 (81) 


which is identical with the expressions of (5.15) (Kielich, 1960). 
Also, eq. (7.5) for the molecular refraction R,, was derived on the simplifying 
assumption of x$ = 1. Without the latter, we res the general expression 


i (i i Ya OF 1 
$^ i dg) as (P + PE i: (32) 
(i) 
which holds for axially symmetric quadrupolar molecules. 


In the case of a one-component system, eqs. (16), (17) and (32) reduce to the 


form 
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4r | 5 9? 
Rm = -z Na h + 3 Gi put 24L pL) + s nis ER ; (35) 
with R denoting the gas constant, 7—the absolute temperature, 6;—the isothermal 
compressibility coefficient of the medium and V—its molar volume. 
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LETTERS TO THE EDITOR 


ON THE HIGH TEMPERATURE SUSCEPTIBILITY OF FERRIMAGNETICS AND 
ANTIFERRIMAGNETICS 


By BOGDAN FECHNER 
Institute of Theoretical Physics, A. Mickiewicz University, Poznań 
(Received October 31, 1960) 


It is well known from statistical mechanics that the zero field susceptibility per atom 


is given by 


kT m 
Zu 2 2 

Kramers has suggested an approximative method of evaluating the function log Z by ex- 
panding it in ascending power series of the reciprocal temperature. The calculations were 
carried through for ferromagnetism and antiferromagnetism by Opechowski (1937) and 
others, for spin 1/2 and by Van Vleck (1937), Brown and Luttinger (1955), Rushbrooke 
and Wood (1958) for arbitrary spin. It is natural to extend these calculations to the case 
of ferrimagnetism and antiferrimagnetism with arbitrary spin of magnetic atoms. In the 
present note we shortly report the results of these calculations. 

We start with the Hamiltonian in the form 


N 
4 -—2 V SSH NS 
D» i 
which relates to an ordered crystal structure composed of two kinds of magnetic atoms. 
Xi, jp means that the summation is extended over all neighbouring lattice sites in which 
magnetic atoms are present. Using the procedure explained by Rushbrooke and Wood, 


we obtain the expansion in the form 


i PPS a 
AT SET LT 


JU 


chere S= n Sat np. CE . The coefficients a, are quantities requiring explicit caleulations, 
= ną: i 
he details of which will he published later. We give here only the general expressions for 
he a coefficients up to n = 3 valid for any crystal structure composed of N4, — A and 
is ' lattice sites occupied by magnetic atoms. 
Vg — B atoms. N4 + Ng = N is the number of lattice sites occupi ; y 8 HP 
Jenoting the three kinds of exchange interactions by J44, Jgg and Jap we fmc 
5 o 


R et 9 S c 1 
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19 = 3 S {na caa [4 (Cag — D Sa — 3] SĄ J34 + np cep [4 (cog — 1) Sp — 3] 58 Jan + 
+ 2n4 Cap [2 (cga — 1) S, + 2 (cag — 1) Sp — 3] 54 Ss Jas + 


ła wę 
+ Bny Caa CAB SA SB Jaa Jap + On 4 CBB CAB DA SB JBB Jap}: 


lą = 135: ee {n4 [2 (Scaa (aa — 1)? — 3e44 — 30p34) Sa — 3(c44 (5c44 — 3) + 10p, 4) Sa + 


+ 10 pa) Sp + 6egg] 53 Jig + na Cap [4 (5 (Cap — 1) (cza — 1) — 3) 6, Sg - 


Es iocus D S 4 ea Cer d S SEIS ud E 
+ 5n4 c44 Cag [A (c44 — 1) 54 — 3] SĄ SB Jaa Jan + 


+ 574 CgB Cap [4 (cgg — 1) Sg — 3] S4 $5 Joe Jap — Na [(6P24,8 + 3C44 aB) + 


t c 90 2 
r 4 (2 Po4,B — C4A CAB (CBA — 1) Sal SA SB Jaa JAB — 
— Ong [(6PaB„A + 3cgB Cpa) + 4 (2 Popa — CBB CBA (Ca — ) Sg] Sa SM Tap ds 
+ 204044 CBB CAB SA SB Jaa J pp Jap}: 


Here 5, stands for 55; (S; + e and n4 = N4/N, ng = N,/N are the concentrations of the 
two components. The p's and c's are structural parameters which have the following meanings: 
744 is the number of nearest A neighbours to an A atom, c45 is the number of nearest B 
neighbours to an 4 atom ete. N 4 ps 4, N 4 ps4 p. NB pag, 4and Ng psp are the numbers of triangles, 
the sides of which join neighbouring lattice sites occupied respectively by: three A atoms, two A 
ind one B atoms, two B and one 4 atoms, three B atoms. The values of these parameters 
for a given structure can be obtained by direct counting on a model lattice. In the special 
case, when there is only one kind of atoms present (Na = N, Ng = 0) our results reduce 
to those obtained by Brown and Luttinger (1955). 


We have also obtained the series expansion of the inverse susceptibility 


ST x 
— E 
Fam Im ETY 


bx 


but the coefficients b, will not be given here with respect to their complexity. They are 
connected with the a coefficients by the relations 


= — — 2 € 3 
by = l, b, = —a,, b, = a4 — ay, b = 20404 — ai — as. 


The application of these results to ferrites and ferromagnetic alloys will be given in 
a subsequent paper. 
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ON A NEW TYPE OF CONDENSER FOR MEASURING THE DIELECTRIC PERMIT- 
TIVITY OF SOLIDS BY THE METHOD OF IMMERSION 


BYL PIECH 
I. Chair of Physies, Mining and Metallurgical Academy at Cracow. 


(Received September 19, 1960; translated paper received on October 21, 1960) 


The “method of immersion” constitutes one of the most exact methods of measuring 
the dielectric permittivity of substances in the solid state. This method, first proposed 
by Starke (1897), is at present applied in all cases of a material accessible in the form of 
minute fragments only, such as powder or small fibres. It bases on the fact that the capacity 
of a condenser filled with liquid does not change if powder of the same dielectric permittivity 
is added to the liquid. Since it is practically impossible to chose the liquid so that it shall 
possess the same permittivity as the powder, the procedure in most cases is this: the con- 
denser is filled in turn with liquids of different values of the permittivity and the powder 
investigated is added, after which the capacity of the system obtained is measured. The 
respective values will be greater or smaller than those for the condenser filled with pure liq- 
uids, according to whether the permittivity of the powder is greater or smaller than that of 
the liquid. By plotting the capacity values of the condenser with the pure liquids on the 
axis of abscissae, and their variations on addition of the powder on the axis of ordinates, 
a graph connecting the points is obtained that intersects the axis of abscissae at the point 
corresponding to the capacity of the condenser filled with a liquid of exactly the same 
permittivity as the powder. With a condenser previously calibrated, the permittivity of 
the powder investigated can be computed. 

The new type of condenser built at this laboratory admits of the rapid and convenient 
determination of the capacity when filled with a pure liquid or with a mixture of liquid and 
powder. Essentially, it presents a system of two cylindrical condensers of approximately 
the same capacity, placed axially one on top of the other (Fig. 1). The metal plates (p) are 
connected to the glass parts (g) with **Araldit", a heat-tempered putty produced by “Ciba”. 
The outer plates are earthed; the inner ones have special leads to the measuring device 
(i, and ij). On removing the screw (s), the powder is introduced and, subsequently, as 
much liquid as is needed for the mixture to fill condenser I and the pure liquid to fill condenser 
II is added. The capacity of condenser II with the pure liquid is now measured. The device 
is now reverted by 180? so that the powder is transferred from condenser I to II, and the 
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capacity of the latter filled with the mixture is measured. Adding a liquid of different permi 
tivity to the condenser, we proceed similarly until the sign of the difference between tl 
capacity of condenser II filled with the liquid and that of condenser I filled with the mixtu 
changes. By the graphical method already described, the capacity is finally measured fi 
a liquid possessing a permittivity value equal to that of the powder. Another procedu 


95 


is the following: if the capacities of the empty condensers I and II are known, these can 
be filled and the the permittivities of the liquids and of their mixtures with the powder 
computed. By plotting the permittivities of the liquids (abscissae) and the differences 
between the permittivities of the respective mixtures and those of the liquids (ordinates), 
a graph is obtained whence the permittivity of the powder is derived directly. 

This type of condenser presents the additional advantage of eliminating error due to 
solubility of the powder in the liquid. Notwithstanding the fact that in such measurements 
liquids are used that do not dissolve the solid investigated noticeably, there is practically 
always some solubility, however small. Hence the use of a single condenser, whose capacity 
is measured once on filling it with the pure liquid and then on introducing the mixture, can 
give rise to error from the change in the dielectric properties of the liquid itself due to 
some amount of the powder in solution. In the present case, however, this can no longer 
occur, as the capacity is measured when the condenser is filled with a liquid that has already 
exerted whatever dissolving properties it possessed. — 

Moreover, a condenser was built for measurements of the the temperature dependence 
of the permittivity (Fig. 2). A system of two condensers remembling the one described is 
surrounded by a double-walled mantle circulating liquid of a temperature stabilized with 
a thermoregulator. Here, the condenser plates are provided by the walls of the hollow 
cylinders, through which leads pass from the internal plates. In order to avoid changes in 
the cápacities of the leads, a system of fixed switches (r) as described by Jezewski (1955) 
is applied. 

'The foregoing types of condensers have been tested at this laboratory by Mr. Mie- 
czysław Wierzbicki and Mrs. Teresa Morstin, who used solutions of nitrobenzene in benzene. 


Starke, H., Wied. Ann:, 60, 629 (1897); 60; 804 (1897). 
Jeżewski, M., Acta phys. Polon., 14, 483 (1955). 
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REVIEWS OF BOOKS 


WiLHELM H. WESTPHAL 


Physikalisches Praktikum, Friedr. Vieweg & Sohn, 
Braunschweig 1959, 285 Seiten, Preis DM. 19,50. 


Physikalisches Praktikum ist eine Sammlung von Beschreibungen praktischer Übungen 
aus dem Gebiete der Physik und ist als Hilfsmittel für die Studenten verschiedener Typen 
der Hochschulen vorgesehen. 

Dieses Buch besteht aus zwei Hauptteilen. Im ersten kurzen Teil (25 Seiten) beschreibt 
der Verfasser die Gróssen, die in den einzelnen Gebieten der Physik auftreten , die Methoden 
der Berechnung der Messungsresultate, gibt eine Reihe praktischer Hinweisungen und Regeln 
betreffs der Messungen. Auch die Fehlerrechnung ist kurz zusammengefasst. Die Aufgabe 
dieses Teiles des Buches ist die Studenten zur korrekter Ausführung und Deutung physikali- 
scher Messungen vorzubereiten. 

Zweiter Teil, beträchtlich linger (250 Seiten), enthält die eigentliche Beschreibung 
von 47 praktischer Übungen. Diese Übungen sind nach folgenden Gebieten der Physik 
eingeteilt: 1) Mechanik 2) Wärme, Gastheorie, Akustik 3) Optik 4) Elektrizität und 
Magnetismus. Am Ende des Buches sind Tabelle für manche physikalische Gróssen in der 
Form der Diagrame angegeben. 

Westphal's Physikalisches Praktikum hat sich besondere Stellung als Handbuch in 
der Praxis im physikalischen Laboratorium erworben. Es ist ein wichtiges Hilfsmittel be- 
sonders aus diesem Grunde, weil der Verfasser sich nicht nur zur Beschreibung der Praxis 
beschrünkt, sondern auch die theoretische Einführung bei jeder Übung angibt. Es ist klar, 
dass die Theorie der Übungen nicht zu ausführlich sein kann, damit das Buch, welches 
den praktischen Charakter haben soll, nicht seine Klarheit verliere. Hier besteht aber 
eine gewisse Gefahr. Im allgemeinen muss die Theorie in einem praktischen Lehrbuch aus 
Notwendigkeit nur skizziert werden und für einen weniger orientierten Leser kann diese 
Tatsache den Eindruck einer allgemeinen Auffassung des Problems hervorrufen. Westphals 
Buch hat eben diese Eigenschaft, dass es einen solchen Anschein nicht gibt, weil der Bereich 
des besprochenen Materials sich sehr eng mit dem Thema der Übung verbindet. Zugleich 
geht aber diese theoretische Einführung in die Einzelheiten des Problems ein und über- 
schreitet oft das Niveau von meisten physikalischen Lehrbüchern für Hochschulen, so dass 
die Beschreibung des Problems nicht verengert wird. Physikalisches Praktikum hat noch 
einen Vorzug: nach der Besprechung jeder Übung folgt ein konkretes Zahlenbeispiel mit 
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den Diagramen und der Diskussion der Fehler, was nicht nur zur besseren Verstindigung 
des Themas beiträgt, sonder auch als Muster bei der Ausführung der Übungen dienen kann. 
Das hilft dem Studenten die fehlerfreie und systematische Arbeit im Laboratorium zu 
erlernen. 

Das Buch hat aber auch einige Nachteile. Der Verfasser hat sich eng auf das Material 
der klassischen Physik beschrünkt. Es handelt sich nicht natürlich um spezielle Übungen 
aus dem Gebiete der modernen Physik, z.B. aus der Kernphysik, aber um gewisse Gebiete, 
welche sich bestimmt in den Bereich der allgemeimen Physik einschliessen lassen. Es wurden 
hier z.B. das Photoeffekt, Optische- und Róntgenspektren und Vakuumtechnik in Betracht 
genommen. Anderseits sind im Physikalischen Praktikum Übungen besprochen, wie die 
Jolly-Waage, welche in der gegenwürtigen Praxis im Laboratorium keine Anwendungen 
finden und welche durch Übungen aus den oben erwühnten aktuellen Gebieten der Physik 
ersetzt werden kónnen. 

Obwohl der Bereich des Materials nich ganz den modernen Bedürfnissen des Laborato- 
riums für Studenten entspricht, so gehórt es doch zu den besten praktischen Lehrbüchern 
und deshalb darf dieses Buch in keinem Laboratorium für Studenten fehlen. Der beste 
Beweis der Vorzüge der guten Bearbeitung des Materials dieses Buches ist die Tatsache dass 
sie seit der Erscheinung (erste Ausgabe im Jahre 1937) schon neun Ausgaben hatte. 

Man soll auch die besonders sorgfiltige graphische Ausführung des Buches hervor- 
zuheben. 


Zofia Leś 
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FRANK MarTosst 
Der Raman-Effekt, Second revised edition. Friedr. Vieweg u. Sohn; Braunschweig 1959, 84 p., 24 figs., 9 tables. 


Matossi' s book was originally edited in the series ,, Verfahrens- und Messkunde der 
Naturwissenschaft". Its primary aim was thus to present experimental and theoretical 
methods employed in Raman spectra investigations. Chapter 1 deals descriptively with 
the fundamentals of the effect, on the basis of the quantum model. In Chapter 2, the author 
considers in turn light sources, Raman spectrographs, methods of measuring the frequencies, 
intensities and degree of depolarisation of the lines, apparatus, and models yielding the 
oscillations of molecules. This last problem initiates Chapter 3 (applications); here, oscillation 
forms for molecules of various symmetry are derived theoretically, and the conclusions 
from Placzek’s polarizability theory are discussed. The oscillations of organic molecules 
and of crystals are treated in some detail. Among the Tables, those of the oscillation types 
and oscillation type characters of molecules of various symmetry are noteworthy. 

The second edition brings much new and valuable material not to be found in the first 
edition published 15 years ago. However, some of this fails to deal adequately with some of 
the most important contemporary problems in Raman techniques. In Chapter 1, a discussion 
of the resonance Raman effect has been added. Chapter 2, although now mentioning the 
loronto type spiral lamp, is totally silent on other types of low-pressure lamps coming 
into use. Some mention is made of photomultiplier methods of spectrum recording; however, 
he author avoids dealing extensively with these growingly employed methods, and refers 
he reader to Brandmueller and Moser’s article (Z. angew. Physik, 8, 35, 142 (1956); 9, 567 
1957)). A description of new types of Raman vessels is added (the multiple reflexion vessel, 
rarious micro-vessels) and the problem of measuring intensities given more place. However, 
only the difference between the maximum and integral intensity is mentioned, whereas no 
mention is made of the attempts at converting the experimental intensities into real intensities 
which are the basis for determining the symmetrical and anisotropic parts of the polarizability 
ensor. At the opening of Chapter 3, the second edition contains a discussion of the 
nethod of determining molecular oscillations by means of a system of oscillating electric 
jreuits; but for this, Chapter 3 is almost unchanged. 

With regard to the small volume, the problems are dealt with very concisely; the reader, 
1owever, will find 143 references to monographs, various collective publications, and to 
he more prominent original papersthat appeared until the end of 1958. The style is uncompli- 
ated, and is easily understood by the reader whose German is not high grade. 
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Matossi's book fails to give an account of the theoretical foundations of the Raman 
effect, whether classical or quantum. This is all the more to be regretted as the applications 
(Chapter 3) are of a predominantly theoretical character (the determination of oscillation 
forms) and are less preoccupied with direct chemical applications. There is not a single 
picture of a Raman spectrum to be found in the entire book. Besides, numerous papers by 
Soviet scientists, especially papers that have not been translated into the Western languages, 


are left unmentioned. 


Roman Mierzecki 


